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Abstract 

A previously  designed  Extended  Kalman  Filter,  based 
upon  the  proportional  guidance  law  and  an  aerodynamic  drag 
equation,  is  modified  to  include  complete  modeling  of  the 
autopilot  and  the  guidance  system  of  the  missile.  A 
sensitivity  analysis  of  three  key  parameters  is  performed 
to  determine  the  observability  of  these  parameters  and  the 
possible  need  to  estimate  them  in  less  complex  filter  models. 
This  analysis  establishes  the  advantages  for  reasonable 
estimates  of  all  three  parameters  during  a high-g  trajec- 
tory. 

A less  ' Dmplex  model  reduces  the  number  of  filter 
states  from  eleven  to  eight.  The  seeker  dynamics,  the  guidance 
system,  and  the  autopilot  of  a missile  are  approximated 
by  a first  order  lag.  The  three  parameters  estimated 
include  the  time  constant  of  the  first  order  lag,  the  pro- 
portional navigation  constant  of  the  missile,  and  the  ratio 
of  mass  to  cross-sectional  area  of  the  missile. 

Under  the  high-g  scenario,  reasonably  accurate 
estimates  of  the  states  and  the  parameters  are  accomplished. 
Under  the  low-g  scenario,  parameter  estimates  are  degraded 
due  to  lack  of  excitation.  This  performance  of  the  para- 
meters, however,  does  not  hinder  the  accurate  estimation 
of  the  missile  states. 
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AN  EXTENDED  KALMAN  FILTER  FIRE  CONTROL  SYSTEM 


! I AGAINST  AIR-TO-AIR  MISSILES 

I.  Introduction  and  Background 

I 

Introduction 

Current  research  into  the  defense  of  an  aircraft  against 
an  air-to-air  missile  has  proposed  the  use  of  an  Extended 
Kalman  Filter  within  the  control  loop  of  an  automatic  fire 
control  system.  An  Extended  Kalman  Filter  (EKF)  in  such  a 
system  would  process  incoming  measurements  of  missile  range, 
range  rate  and  line-of-sight  angle  to  achieve  both:  current 
state  estimation  - to  aid  in  pointing  accuracy,  and  for 

f 

j trajectory  prediction  - to  aid  in  probability-of-kill 

assessment.  The  capability  of  the  Extended  Kalman  Filter  in 
this  role  is  directly  dependent  upon  the  accurate  modeling  of 
the  dynamics  of  the  missile.  The  more  accurate  the  model 
of  the  missile,  the  better  the  estimation  and  prediction 
capability  of  the  Extended  Kalman  Filter. 

Early  attempts  at  state  estimation  of  air  targets  used 
a first  order  Gauss-Markov  process  as  an  acceleration  model. 
An  appropriate  choice  of  process  noise  strength  and  correla- 
tion time,  along  the  three  axes  of  movement,  was  used  in 
modeling  the  target's  acceleration  components.  The  results 
of  these  studies  have  consistently  implied  that  the  filter's 
estimates  were  very  sensitive  to  the  least  mismodeling  of  the 
target's  acceleration.  The  model  is  generic  in  nature  and 
has  provided  satisfactory  results  in  several  cases 
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(Ref  2:366;  8:323).  However  to  apply  this  approach  in  model- 
ing an  air-to-air  missile,  ignores  two  assumptions  which  can 
be  made  about  this  type  of  air  target: 

1)  most  air-to-air  missiles  use  proportional 
navigation  as  a guidance  system  and, 

2)  many  current  air-to-air  missiles  are  "boosted 
to  Mach"  and  coast  to  intercept. 

In  current  thesis  work  at  the  Air  Force  Institute  of 
Technology,  an  alternate  approach  has  met  with  moderate 
success  by  utilizing  these  two  assumptions  (Ref  5) . These 
results  were  incomplete,  but  definitely  established  the 
feasibility  of  this  approach  and  demonstrated  the  advantages 
of  incorporating  these  two  assumptions  into  a model  of  the 
missile's  acceleration.  Essentially,  the  model  attempts  to 
exploit  more  physical  knowledge  about  the  missile's  accelera- 
tion than  the  generic  Gauss-Markov  model.  This  study  extends 
the  results  obtained  by  Lutter  in  Reference  5 by  also  using 
these  two  assumptions. 

Background 

Proportional  navigation  has  proven  to  be  straight- 
forward to  implement  and  tactically  effective.  Although 
this  guidance  law  is  not  optimal  for  most  realistic 
intercept  problems,  those  attempts  to  implement  optimal 
guidance  on-line,  usually  approximate  proportional  navigation 
in  some  form  (Ref  12) . Also,  because  of  the  inherent 
simplicity  to  mechanize  proportional  navigation,  it  seems 
likely  that  this  guidance  law  will  be  retained  for  some  time. 
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particularly  in  air-to-air  missiles  (Ref  11:18). 

A guidance  system  with  proportional  navigation  attempts 
to  null  the  angular  rate,  9,  of  the  line-of-sight  (LOS)  vector 
between  the  target  (ownship)  and  the  interceptor  (missile) . 
This  is  performed  by  the  missile  commanding  a lateral 
acceleration  proportional  to  closing  velocity  and  the 
amgular  rate  of  LOS.  Analytically, 

= n''J  a) 

where 

a = commanded  missile  lateral  acceleration 
c 

n = proportional  navigation  constant 

V = closing  velocity  along  the  LOS  vector 

w 

0 = angular  rate  of  LOS  vector  with  respect  to 

inertial  space. 

The  second  hypothesis  of  a non-thrusting  missile  implies 
that  acceleration  along  the  body  axis  of  the  missile  can  be 
well  approximated  by  the  aerodynamic  equation 

where 

a^j  = the  missile  acceleration  due  to  drag 
(deceleration) 
p = the  air  density 
Cjj  = the  coefficient  of  drag 

V * the  velocity  of  the  missile 
m 

S = the  missile  cross-sectional  area 


M = the  mass  of  the  missile. 


As  with  most  engineering  problems,  the  benefits  gained 
by  a particular  approach  are  usually  tempered  by  some  added 
complexity.  In  this  case,  the  use  of  Equations  (1)  and  (2) 
require  knowledge  of  three  parameters;  C^,  n,  and  M/S  for 
implementation.  An  expression  for  Cj^  in  terms  of  measure- 
able  quantities  was  developed  by  Lutter  (Reference  5)  which 
well  approximates  the  coefficient  of  drag  for  a generic 
missile  at  speeds  above  Mach  1.1,  This  approximation  is  not 
very  restrictive  since  the  velocities  of  air-to-air  missiles 
are  well  above  this  speed  in  all  combat  scenarios.  The 
parameters  n and  M/S  may  either  be  chosen  using  a priori 
knowledge  of  an  attacking  missile  or  may  be  estimated  by 
using  the  unknown  parameter  estimation  capability  of  an 
EKP.  The  former  approach  includes  the  use  of  nominal  values 
for  these  parameters.  The  latter  choice  is  accomplished 
by  augmenting  the  state  model  of  the  missile  in  the  filter 
with  additional  states  modeling  the  unknown  parameters.  The 
estimation  of  the  two  parameters  n and  M/S,  although  adding 
complexity  to  the  problem,  can  provide  useful  information 
for  probability-of-kill  assessment  and  threat  identification 
to  the  fire  control  system.  Lutter  chose  to  estimate  n and 
M/S. 

To  conduct  a meaningful  performance  analysis,  Lutter 
developed  a realistic  model  of  a missile.  This  "truth 
model"  included  significant  detail  of  the  missile's  kine- 
matics, dynamics,  and  internal  structure.  Specifically, 
the  guidance  strip,  shown  in  Figure  1,  modeled  the  lag  and 


filtering  effects  of  the  three  components  shown.  However, 
in  order  to  reduce  the  complexity  of  the  model  in  the  filter, 
Lutter  ignored  these  effects  and  equated  commanded  accelera- 
tion (defined  in  Equation  (1))  to  developed  lateral 
acceleration.  The  results  of  his  study  indicated  the  need 
for  a more  refined  model  of  the  guidance  strip. 

Statement  of  the  Problem 

The  problems  addressed  by  this  thesis  work  are  the 
following ; 

(1)  The  feasibility  of  the  fundamental  Filter 
design  has  already  been  established,  but  a benchmark, 
showing  the  best  estimation  possible,  was  not  considered  in 
previous  work.  The  first  topic  addressed  will  determine  a 
proper  benchmark  which  will  be  used  for  comparison  with 
"reduced  order"  filters  ("reduced  order"  referring  to  the 
number  of  states,  of  the  guidance  strip)  . The  designing  of 
the  benchmark  will  emphasize  the  key  components  of  the 
guidance  strip. 

(2)  A sensitivity  analysis  will  be  performed  for 
those  parameters  which  are  not  normally  known  a priori,  or 
cannot  be  defined  in  terms  of  observable  quantities,  and 
must  be  known  for  reduced  order  filters.  An  excimple  of 
this  is  the  proportional  navigation  constant,  n.  This 
analysis  will  better  define  the  tradeoff  between  a nominal 
value  or  estimating  each  unknown  parameter. 

(3)  A reduced  order  filter  will  be  designed  and 
evaluated.  To  model  the  complete  guidance  strip,  assuming 
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that  all  parameters  are  kncwn,  would  require  five  filter 
states.  The  approach  taken  by  Lutter  was  to  model  this  with 
a zero  order  transfer  fiinction.  The  insights  gained  from  the 
benchmark  design  and  the  sensitivity  analysis  will  influence 
proposed  reduced  order  designs  of  greater  complexity  than 
the  zero  order  model.  The  minimum  design  to  be  considered 
will  be  a first  order  lag. 

Monte  Carlo  Analysis.  The  accurate  description  of  the 
missile  kinematics  in  a proper  filter  model  required  the 
use  of  nonlinear  differential  equations.  Therefore,  the 
conventional  Kalman  Filter  was  not  an  appropriate  estimator 
for  this  particular  system.  The  Extended  Kalman  Filter 
(EKF) , as  chosen  by  Lutter,  was  a good  initial  choice  as  a 
means  of  solving  this  nonlinear  estimation  problem.  The 
evaluation  of  the  EKF  requires  the  use  of  a Monte  Carlo 
analysis.  A complete  description  of  a Monte  Carlo  analysis 
is  included  in  Reference  6.  In  general  this  analysis 
involves  an  ensemble  average  of  the  filter  error  over 
a specified  number  of  simulation  runs.  The  nimnber  of 
runs  chosen  is  usually  dictated  by  the  confidence  required 
in  these  compiled  statistics.  A related  problem  in  this 
thesis  work  will  be  choosing  the  approriate  number  of  runs 
to  obtain  results  with  reasonable  confidence. 

Software  Changes.  The  software  provided  by  Lutter  was 
modified  to  suit  the  needs  of  this  particular  study.  The 
modifications  and  reasons  for  each  are  included  in  Appendix 
A.  The  results  of  the  Lutter  design  appearing  in  this 
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thesis  include  those  changes. 
Organi zation 


The  body  of  this  thasis  las  been  arranged  into  four 
additional  chapters.  Chapter  suvjaarizes  the  key  ideas 

and  equations  used  by  Latter  in  developing  tiie  truth  ncdel 
for  this  problesi.  Ch.apter  Three  includes  the  development 
of  the  basic  filter  derived  by  Lutter.  In  addition , the 
benchmark  and  reduced  order  filters  .gerenatod  in  this  study 
are  preronted.  The  sensitivity  aralysis  of  the  unknown 
parameters  as  well  as  the  estimation  ani-lysis  of  the  reduced 
order  filter  with  parameter  estimation  comprise  Chapter 
Four.  The  final  results,  conclusions  and  suggestions  for 
further  analysis  conclude  this  \;ork  in  Chaptei  Five. 

Appendix  B was  constructed  to  contain  ail  plots 
produced  during  the  analysis.  It  was  necessary  to  make  this 
a separate  volume  becuase  of  the  large  :._mb;er  of  plots 
generated  from  the  myriad  of  changes  performed  during  the 
filter  analysis.  However,  the  important  plots  are  contained 
within  the  body  of  this  report  vrhen  necessary  as  supportive 
material . 
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II.  TRUTH  MODEL 


Introduction 

Oversimplification  in  modeling  the  "true  dynamics"  of  a 
non-linear  system  can  degrade  the  analysis  of  an  Extended 
Kalman  Filter  and  misrepresent  its  actual  performance  capa- 
bilities. The  attention  to  detail  should  be  an  overall 
qualifier  for  the  results  of  the  simulation.  With  this 
thought  in  mind,  the  "truth  model"  (those  equations  defining 
the  system  of  interest)  has  been  modeled  to  include  signifi- 
cant detail  in  the  simulation  of  the  kinematics,  dynamics, 
and  internal  structure  of  an  air-to-air  missile.  The  truth 
model  also  includes  a flexible  model  of  an  ownship  target 
which  simulates  the  maneuvers  of  an  evading  aircraft. 

The  following  chapter  presents  the  coordinate  frames  and 
the  engagement  geometry  of  the  truth  model  developed  by 
Lutter  (Ref  5:7-50).  Also  presented  are  the  most  important 
equations  of  the  truth  model  with  the  basic  assumptions  used 
in  their  development.  However,  the  presentation  in  this 
report  is  brief  and  serves  only  as  a review  of  the  develop- 
ment by  Lutter.  Those  readers  desiring  a complete  develop- 
ment should  refer  to  this  earlier  development  (Ref  5:7-50). 
Coordinate  Frames 

The  filter  was  designed  with  the  intent  of  applying  it 
to  cases  which  involve  multiple  targets.  This  design  goal 
led  to  the  use  of  a local  inertial  frame  which  reduces  the 
computational  burden  that  would  be  necessary  with  the  use 
of  multiple  line-of-sight  frames  (Ref  5:1).  The  analysis 


is  restricted  to  a two  dimensional  scenario  in  the  horizontal 
plane.  The  local  inertial  frame  is  defined  to  lie  parallel 
to  this  horizontal  plane  of  motion.  The  origin  is  arbitrarily 
chosen  to  be  the  position  of  the  ownship  at  the  start  of  the 
intercept . 

In  addition  to  the  inertial  frame,  the  geometric  develop- 
ment exploits  three  other  coordinate  frames.  These  coordi- 
nate frames  are 

1)  the  ownship  tracker  line-of-sight  frame  (TLOS) , 
used  in  conjunction  with  the  measurements  of  the  missile 

2)  the  initial  line-of-sight  frame  of  the  missile 
(IMLOS) , used  to  define  the  missile  kinematics  within 
the  truth  model  and  to  measure  the  angular  rates  of  the 
line-of-sight  vector  to  the  ownship  (maneuvering 
aircraft) 

3)  the  initial  tracker  line-of-sight  frame  (ITLOS)  used 

to  relate  missile  motion  in  the  IMLOS  frame  to  the  local 
inertial  frame,  defined  such  that  "^IMLOS,  ^ITLOS 

= “^IMLOS. 

Figure  2 shows  the  relationship  for  these  four  frames  at 
time  (t)  = 0 and  t = t^ . Note  that  at  t = 0 the  ITLOS 
and  TLOS  are  aligned.  Also,  at  t = 0,  the  missile  line-of- 
signt  angle,  0,  is  defined  as  zero. 

Geometry 

The  engagement  geometry  may  at  first  appear  complicated 
because  of  the  many  angles  defined  but,  once  understood,  it 
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can  be  appreciated  for  its  direct  approach  and  versatile 
nature.  Figure  3a  describes  the  complete  engagement  geometry 
at  the  beginning  of  the  intercept  (t=0) . This  figure  relates 
the  intercept  geometry  of  the  missile  to  the  engagement 
geometry  of  the  truth  model.  In  this  figure  the  initial 
missile  velocity  and  initial  ownship  velocity  are  represented 
by  the  vectors  and  V^,  respectively.  The  missile  lead 
angle,  6^^,  defines  the  angle  between  the  missile's  current 
velocity  vector  and  the  X-axis  of  the  IMLOS  frame. 

In  order  to  actively  estimate  the  change  in  the  LOS 
vector  that  the  missile  observes,  the  angle  6^  is  defined 
(Figure  3a,  Figure  3b) . Using  a basic  law  of  geometry, 
that  opposite  interior  angles  formed  by  a line  intersecting 
two  parallel  lines  are  equal,  the  following  equation  is 
formed  (Figure  4) 

0^  - 180  = 0 + ().  (3) 

where 

0^  and  0 are  as  previously  defined 

(fr  = the  angle  through  which  IMLOS  frame 
is  rotated  from  the  local  inertial 
frame  (constant  in  time) 

Differentiating  Equation  (3)  with  respect  to  time,  the 
desired  relationship  for  the  angular  rate  of  the  LOS  vector 
becomes  apparent. 

0^  = 6 (4) 
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(original  position) 


Thus,  variations  in  9^/  measurable  at  the  ownship,  reflect 
those  variations  seen  by  the  missile  in  0 . Recalling  the 
model  of  the  proportional  guidance  law  to  be  used  in  the 
filter,  the  merits  of  measuring  6^,  in  this  manner,  becomes 
quite  obvious; 

/ 

a^  = nv^e  (1) 

This  angle,  0^,  along  with  range  and  range  rate  will  be  the 
three  measurements  passed  to  ti.e  Extended  Kalman  Filter  for 
processing. 

Figure  5 depicts  the  typical  intercept  geometry  for  an 
air-to-air  missile.  The  angular  rate  of  the  LOS  vector  to  the 
ownship  is  measured  by  the  missile  with  respect  to  some 
arbitrary  non-rotating  reference  for  missiles  with  inertially 
stabilized  seekers.  In  this  case,  the  reference  chosen  is 
the  LOS  vector  that  is  first  sighted  by  the  missile  upon 
beginning  intercept.  The  origin  of  the  IMLOS  frame  is 
defined  to  be  the  center  of  the  seeker  antenna  of  the  missile. 
The  lateral  acceleration,  aj^,  of  the  missile  is  defined  as 
positive  when  causing  an  increase  in  the  lead  angle  of  the 
missile,  d^.  This  relationship  is  also  true  for  the 
acceleration  of  the  ownship,  a^,  and  the  angle,  0^.  0^  is 

defined  as  the  angle  between  the  axis  and  the  ownship 

velocity  vector.  Missile  guidance  controllers  are  usually 
implemented  with  a small  angle  approximation  for  the  LOS 
angle,  0,  where  9 is  well  approximated  by  the  ratio  of  the 
lateral  distance,  y^,  to  the  range,  R.  The  missile  truth 
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Fig.  5.  Intercept  Geometry  for  Missile 
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model  in  this  simulation  does  include  the  sn’ell  angle  approyi- 
m.ation  for  0 


6 = arcsin  (y^/R)  = 


The  intercept  geometry  for  an  air-to-air  missile  is  used 
to  define  for.  stat'-'  variables  of  the  truth  riiodelf  all  quanti- 
ties being  referenced  with  respect  to  the  IMLOS  frame. 


_ V„sin0,  - V sinO^ 
^ d = T ii  m L 


A = a^,g/V,^, 


' - ] = -2 

K -V  R 

c 


where 


0 = the  angle  between  the  axis  and  the  ownship 

velocity  vector 

0_  = the  angle  between  the  X axis  and  the  missile 

JLi 

velocity  vector 

= the  magnitude  of  the  velocity  vector  of  the  ownship 

V - tlK-  maqnitudc  of  the  v'elocity  vector  of  the  missile 
= the  tracker  .acceleration  perpendicular  to 

expressed  in  g's 

V = the  realtive  closing  velocity  between  the  missile 


and  the  trac}:or 


R and  as  defined  previously 


IHIS  IS  .£0  DOfi 
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The  orig5.n  of  the  ownship  103  frarie  (TLOS ) is  defired  r 
the  center  of  its  inea-.uring  antenna.  The  antenna  of  Che  ov/n- 
chip  is  assumed  to  point  withouc  error  and  v/ithout  lag  effects. 
This  simplification  assumes  that  the  possible  errors  arising 
from  the  dynamics  of  the  antenna  are  negligible  as  compared 
to  the  measurement  u certainty. 

Thus  far,  three  coordinate  frames  have  been  discussed: 

1)  Local  inertial  frame 

2)  IMI.03 

3)  TLOS 

Since  the  filter  estimates  v;ill  be  described  with  respect  to 
the  inertial  frame,  the  true  missile  velocity  and  position 
must  also  bo  defined  with  respect  to  the  inertial  frame  for 
an  estimate  of  true  error.  This  is  accomplished  by  relating 
the  current  missile  acceleration  vectors  in  the  IMLOS  frame 
to  the  ITLOS  frame  (recalling  that  = _^IML0S  ,^ITL0S  ^ 

, IMLOS)  : 


^ITLOS  ^ a cos6  + a..  sinG 
mx  D L 1 L 


^IPLOS  _ 3 _ a cos 9 

my  D L L L 


v;herG 


= the  missile  drag  acceleration  (deceleration) 

= the  lateral  acceleration  of  the  missile 

Li 

„ITLOS  ,, ITLOS  1-4.  X,  • ■ 1 • 4.V, 

V , , V the  velocitv  of  the  missile  in  the 

iTiX  my 

ITI.OS  frame  in  the  x a)id  y direction 


rrspectivc')  y 


These  quantities  are  then  ooordinitized  into  the  local  inertial 
frame  using  the  initial  9 defined  as  0^. 


v^3^°s^ose. 

mx  0 

my 

^ITLOSg^  Q 

..ITLOS^^^e 
— V cos 

mx  0 

my 

0 

0 


(12) 

(13) 


where 

and  = the  velocity  of  the  missile  in  the 
mx  my 

inertial  frame  in  the  x and  y direction 
respectively. 


Truth  Model 

Figure  6 represents,  in  block  diagram  form,  the  overall 
truth  model.  Each  block  represents  the  logical  grouping 
of  equations  modeling  a specific  function  which  is  addressed 
sequentially  as  they  appear  in  the  figure. 

Ownship  Kinematics.  The  motion  of  the  ownship  is 
described  by  the  following  equations : 


kj^t  + k2COs  (ut) 

(14) 

kjt  + k^sin((Dt) 

(15) 

ki  - k2a)sin(a)t) 

(16) 

k^  + k^a)cos(a)t) 

(17) 

2 

(18) 

2 

-k^co  sin  (wt) 

(19) 
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Pig.  6.  Block  Diagram  of  the  Truth  MacJel 
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where 


k2»  kj,  k^  (^  = varic±)les  chosen  for  desired 

performance  of  ownship 

These  equations  were  chosen  to  represent  the  motion  of  the 
ownship  because  they  embody  the  characteristic  high-g  turning 
and  jinking  maneuvers  of  an  evading  aircraft.  The  variables 
of  the  equations  were  chosen  more  for  the  resultant  effect 
upon  the  trajectory  of  the  missile,  than  to  model  a particu- 
lar evasion  maneuver.  The  two  sets  of  variables  used  in  the 


analysis  are  shown  in  Table  I. 


Table  I 

Ownship  Variables 


The  first  set  produced  a rather  benign  trajectory  with  the 
missile's  lateral  acceleration  never  exceeding  six  g’s.  The 
second  set  subjected  the  missile  to  a high-g  trajectory  with 
the  missile  reaching  lateral  accelerations  as  high  as  twenty- 
three  g's.  Figure  7 depicts  the  lateral  acceleration  profiles 
for  both  the  low-g  and  the  high  g scenarios.  Figures  8a 
and  8b  present  the  flight  profiles  of  both  the  aircraft  and 
the  missile  in  low  and  high-g  trajectories.  The  marking 
identifiers  on  the  profile  plots  are  positioned  at  constant 
intervals  of  0.56  seconds. 

Missile  Seeker  Noises . This  section  models  the 
dominant  aspects  of  seeker  noise  for  a wide  class  of  air-to- 
air  missiles.  The  uncertainty  corrupting  the  measurements 
of  the  missile's  seeker  can  be  placed  into  three  broad 
categories:  glint,  amplitude  scintillation  and  thermal 

noise.  Those  effects  which  corrupt  seeker  angle  measure- 
ments are  collectively  referred  to  as  glint  (or  angle  scintil- 
lation)'. The  amplitude  variations  of  the  received  signal  are 
aptly  termed  amplitude  scintillation.  Background  noises  and 
system  inaccuracies  are  grouped  together  as  thermal  noise. 

Glint  noise  is  caused  by  the  wandering  of  the 
apparent  centroid  of  radiation  of  the  target  (Ref  13:97). 

These  effects  are  time-correlated  and  are  interpreted  by  the 
missile  guidance  system  as  actual  changes  in  the  angular 
velocity  of  the  target  (ownship) . The  dominant  charac- 
teristic of  a time-correlated  process  such  as  angle 
scintillation  noise  can  be  captured  by  modeling  it  as  a 
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first  order  Gauss-Markov  process,  generated  by  passing  a 
white  Gaussian  noise  through  a first  order  lag  (Ref  6:192). 

The  lag  time  constant  and  input  strength  of  the  white  noise 
must  be  appropriately  chosen  for  the  target  and  system 
\incer taint ies  being  considered. 

Amplitude  scintillation  is  a function  of  target 
aspect  angle.  More  simply,  the  changes  in  the  effective 
cross-sectional  area  of  the  ownship,  facing  the  missile  seeker, 
cause  variations  in  the  amplitude  of  the  received  radiation. 

The  inaccuracies  in  missile  guidance  caused  by  amplitude 
scintillation  are  considered  negligible  for  the  ranges  con- 
sidered in  this  study,  and  for  this  reason  are  not  modeled. 

The  modeling  of  thermal  noise  attempts  to  include 
those  effects  which  are  uncorrelated  in  time  and  are  rela- 
tively wideband  with  respect  to  the  "guidance  strip"  of  the 
missile.  The  "guidance  strip"  of  the  missile  includes 
those  components  used  in  processing  line-of-sight  informa- 
tion,  6' , into  missile  lateral  acceleration.  For  a radar 
seeker,  a good  example  of  thermal  noise  would  be  unexplain- 
adjle  voltage  changes  in  electronic  equipment,  while  for  an 
infrared  seeker,  the  temperature  of  the  background  atmos- 
phere or  under  cover  are  possible  thermal  noise  sources. 

The  thermal  noises  have  been  modeled  as  white  Gaussian  noises. 

In  general,  the  effects  of  glint  and  thermal  noise 
are  range  dependent:  thermal  noise  effects  decrease  and 
glint  effects  increase  as  range  decreases.  However,  for 
this  study,  the  range  dependency  has  not  been  modeled.  A 
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realistic  strength  of  these  noises  for  a radar  system  at  a 
range  of  2,000  feet  is  chosen  and  held  constant  over  the 
entire  intercept  scenario.  Table  II  shows  strengths  of 
the  uncorrelated  and  correlated  noises  as  well  as  the  lag 
time  constant  of  the  correlated  noise. 


Table  II 


Sigmas  and  Time  Constants  for  Seeker  Noise  (Q'-Q) 


Types  of  Noises 


Sigmas  (a) 


tincorrelated  (Thermal)  0.003 


correlated  (Glint) 


0.00894 


Time  Constant 


0.1  (Sec) 


Missile  Seeker.  In  this  analysis  the  seeker  has 
been  modeled  to  perform  the  following  functions: 

1)  track  the  ownship  continuously 

2)  measure  the  LOS  angle  rate . 

A characteristic  of  all  seekers  is  the  distortion  effect, 
aberration  error,  of  the  radome  (or  irdome)  upon  the 
received  radiation.  Figure  9 shows  the  geometry  used  in 
modeling  the  seeker  and  explicitly  defines  the  error,  9^, 
caused  by  aberration  error.  With  aberration  error,  the 
seeker  now  attempts  to  align  the  centerline  of  the  antenna 
with  the  LOS  to  the  ownship  by  nulling  the  false  borosight 
error  t' . Since  the  aberrration  angle  error  varies  with  look 
angle  / there  is  a coupling  of  missile  pitch  rate  with 

measured  ownship  motion.  Figure  10  depicts  a plot  of 
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Fig.  10.  Aberration  Angle  Error  vs.  Look  Angle 
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eUserration  angle  error  as  a function  of  look  angle  (Ref  10;A-6). 
The  dominant  aspects  of  aberration  error  have  been  modeled 


j 

I 


I 


I 
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as  a constant  slope,  of  .001  and  with  an  angle  bias  of 
zero. 

The  dynamics  of  the  seeker  have  been  modeled  as  a 

first  order  lag  with  a time  constant,  of  .075  seconds. 

A block  diagram  showing  the  input  9 and  the  coupling  effect 

that  k^  has  upon  missile  body  motion  with  ownship  motion  is 

shown  in  Figure  11.  This  effect  can  be  seen  by  examining 

the  summing  junction  prior  to  and  noting  the  contribution 

of  0_  to  0 is  canceled  if  k =0.  As  shown  in  this  figure, 

the  output  of  the  seeker  and  track  loop  is  the  measured 

LOS  rate , 6 " . 

s 

Missile  Guidance.  The  output  of  the  seeker  is  a 
noisy  measurement  of  the  LOS  rate.  This  is  true  for  most 
missiles  and  therefore  implementation  of  some  form  of 
filtering  is  necessary  to  reduce  the  noise  on  the  9'  signal. 
Usually  the  filtering  is  performed  by  the  guidance  system, 
with  a variable  gain,  to  produce  a desired  commanded 
acceleration,  a (Ref  11:1).  The  filtering  of  9'  has 
been  modeled  as  a first  order  lag  with  a time  constant 

(T2)  of  .3  seconds.  The  complete  transfer  function  for 

• 

9"  to  commanded  acceleration  is 
s 

V 9'’ 

- ri  c s (20) 

st2  +1 


where  all  variables  have  previously  been  defined.  Although 

31 


somewhat  simplified  in  approach,  this  modeling  does  approxi- 
mately characterize  some  guidance  systems  found  in  modern 
air-to-air  missiles  (Ref  11:42).  The  missile  guidance  system 
will  later  be  shown  to  be  the  dominant  block  of  the  "guidance 
strip." 

Airframe/Autopilot . The  missile  airframe  modeled  is 
a symetrical  cruciform  with  four  fixed  wings  and  four  rear 
moveable  control  surfaces  (Figure  12) . The  cruciform  con- 
figuration permits  lateral  maneuvering  in  any  direction 
without  rolling.  Because  of  this  symmetry,  the  horizontal 
and  vertical  planes  can  be  considered  uncoupled.  The  major 
assumptions  used  in  modeling  the  dynamics  of  the  missile 
are 

1)  the  mass  of  the  missile  is  constant 

2)  the  missile  is  roll  stabilized 

3)  the  missile  is  a rigid  body  and  can  be  represented 
as  a point  mass 

4)  small  angle  approximations  are  valid  for  the 
angle  of  attack. 

The  first  ass\imption  is  consistent  with  the  intent  of  the 
problem  to  consider  the  aerodynamic  forces  acting  upon  the 
missile  after  it  has  been  boosted  to  Mach  and  its  fuel  is 
spent.  The  second  assumption  is  nonrestrictive  since  a wide 
class  of  air-to-air  missiles  are  roll  stabilized  to  reduce  the 
complexity  of  the  guidance  system  (Ref  1:235;  11:23).  The 
third  assumption  considers  all  forces  to  act  through  the 
center  of  gravity  as  well  as  neglecting  all  rotational 
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dyneutiics.  The  fourth  assumption  ignores  the  nominal  values 
of  drag  caused  by  increasing  the  effective  cross  sectional 
area  of  the  missile  into  the  effective  wind.  This  is  true 
because  the  small  angle  approximation  assumes  that  the  effec- 
tive wind  is  always  aligned  with  the  longitudinal  axis  of 
the  missile.  It  also  enables  one  to  make  the  assumption  that 
the  developed  lateral  acceleration  is  perpendicular  to  the 
velocity  vector. 

Given  the  above  assumptions,  the  four  equations 
which  describe  the  dynamics  of  the  missile  are 


= - \ {“(t)  - q(t)} 

(21) 

q(t) 

= M q(t)  + M a(t)  + M 5 (t) 
q ct  5 

(22) 

a(t) 

= q(t)  - L^a(t)  - L^S(t) 

(23) 

Mt) 

= -X6(t)  + Xu(t) 

(24) 

where 

aj^(t)  = missile  lateral  acceleration  due  to  lift 

V_  = missile  velocity  maanitude 

m 

a(t)  = missile  angle  of  attack 

q(t}  = missile  pitch  rate,  note;  0^^^  = q(t) 

5 (t)  = control  surface  deflection 

u(t)  = control  command  (commanded  value  of  6) 

-X  = actuator  pole  (bandwidth) 

M , , M. , L , L . = Missile  stability  derivatives 

q ao  a o ■* 

(Ref  4:14;  5:20-21;  9;  Chap.  8,  p.3) 
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The  control  coroinand  , u(t),  is  r.'.odelcd  as  a function  of 
coitmanded  acceleration,  pitch  rate  and  lateral  acceleration. 


u(t)  = kga^(t)  - h2a^(t)  - hj^q{t) 


(25) 


After  so:ae  algebraic  inanipulation , equations  (21)  through 
(25)  can  be  formed  into  the  system  of  equations 


a (t) 

- (L  th-XV  Ls ) 
a 2 m 6 

(V  T,  -h.XV  L.)  (XV  L,) 

m a 1 mo-  mo 

aj^(t) 

JLl 

g g 

q (t) 

* 

M 

q(t) 

ra  a 

q 

a 

6 (t) 



-gh-,X 

-h^X  -X 

6 (t) 

XV’  lH 
0 m 


gXk, 


left)  (0,5, 


where  the  gravity  unit  g = 32.2  ft/sec  has  been  ineJudod 
for  proper  "cnling. 

The  missile  has  been  modeled  with  an  adaptive  autopi- 
lot which  varies  the  gain  k^  and  feedback  gains  h^,  112  to 
maintain  a specif  j.c  pole  placement  of  the  airframe  response  . 
The  desired  response  is  de.scribed  by  the  second  order  system 


w 2 

_ n 

2Co  S t co^ 
n n 


(27) 
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where 


C = .707 


quality  prachcable 

I*35UW  COi- Y iuiu'ilohJuJ  IV  DUG  __ 


Oi  = 7.07  rad/'s ec. 
n ' 

To  accomplish  this  k - , h^^  and  h^  are  varied  according  to  the 
following  equations 

h,  1 M ^ L M 

. 1 . V,  1 a+  Ot  CT  , r 

kr  — — + h — g — . = — r, \ ^ * 

6 V 2 V L y..  - ^J  L. 

m mao  a o 


, - tI  h,  + '-'a'  <l7  * -V^ 

’’l  ■■  M,  "2 

' i:  L.  + “ ) + (M„  + 2t»)  ILsH„-  H,> 


0 Ci  0 Ci 


h\  «'  ■ - 2;^  - M^)  4-  - q)  - ^ - I-i„) 


The  stability  coefficients  and  stability  derivatives 
are  determined  within  the  simulation  by  use  of  a cubic  fit, 
obtained  from  the  Air  Force  Avionic  Laboratory  (AFAL) , to  a 
set  of  empirical  data  from  actual  flight  tests  of  an  air- 
to-air  missile.  This  method  provides  accurate  estimates  of 
these  coefficients  and  derivatives  for  the  adaptive  gains 
used  in  the  autopilot.  Equation  (26)  vrill  later  be  used 
to  define  the  airframe/autopilot  transfer  function.  Ey 
assuming  time  invariant  coefficients,  the  Laplace  tra.nsform 
of  this  transfer  function  may  be  found  to  be 


^2 


+ ^1^  + ^0 


+ b_j^S  + bQ 
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where 


32#  a^,  aQ,b2,  b^,  are  defined  as  the  airframe/ 
autopilot  coefficients 

During  the  analysis,  these  coefficients  were  calculated  at 
specific  times  to  generate  a model  of  this  transfer  function. 
This  is  discussed  further  in  Chapter  III. 

Missile  Kinematics . The  missile  kinematics  can  be 
described  by  the  two  equations 

" I P I Vm^ 

and 

®L  ' 


(32) 
( 7) 


where  a^j  = acceleration  magnitude  due  to  drag 
p = density  of  the  air 
S = cross  sectional  reference  area 
M = mass  of  the  missile 
Cp  = coefficient  drag 

= velocity  magnitude  of  the  missile 
a_  = lateral  acceleration 

Xj 

g =32.2  ft/sec^ 


Equation  (7)  was  originally  introduced  with  the  interecept 
geometry.  The  coefficient  of  drag  can  be  well  approximated 
for  speeds  greater  than  1.1  Mach  by  the  following: 


Co  - V 

(V_)‘’  Sa 


in 


(33) 
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where 


A = speed  of  sound  at  appropriate  altitude,  air 
temperature 

= coefficient  of  normal  force  with  respect  to 
the  angle  of  attack 
Ct  = coefficient  of  lift  (Ref  5:7-50) 


The  coefficient  of  lift,  , may  be  defined  as 

Xj 

C-  = . M 

^ — 1 S 


(34) 


pv, 


m 


By  substituting  Equation  (32)  into  Equation  (33)  Cj^  is 
expressible  as 


2(A)^  + 

(V„)^ 

m No  r m 


M 

S 


(35) 


The  reference  area,  S,  is  usually  defined  as  the  cross 
sectional  area  of  the  missile.  This  makes  it  easy  to  scale 
for  different  missiles.  From  Equation  (35)  it  is  seen  that 
Cjj  is  a function  of  S.  From  this  realtionship  any  S will 
yield  an  appropriate  when  solving  for  the  a^  in 
Equation  (32).  With  this,  both  equations  of  the  missile 
kinematics  are  completely  defined. 

Tracker  Noises.  The  tracker  being  considered  in 
this  study  is  modeled  as  a typical  radar  for  a fighter  type 
aircraft.  A modeling  of  the  uncertainty  of  the  radar 
measurements  follows  a similar  development  as  was  described 
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where 


= speed  of  sound  at  appropriate  altitude,  air 
temperature 

= coefficient  of  normal  force  with  respect  to 
the  angle  of  attack 
= coefficient  of  lift  (Ref  5:7-50) 


The  coefficient  of  lift,  C^.  , may  be  defined  as 


C_  = . M 

^ s 

m 


By  substituting  Equation  (32)  into  Equation  (33)  is 


expressible  as 


= 2(A)’*  + 

'D  rr  
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The  reference  area,  S,  is  usually  defined  as  the  cross 
sectional  area  of  the  missile.  This  makes  it  easy  to  scale 
for  different  missiles.  From  Equation  (35)  it  is  seen  that 
is  a function  of  S.  From  this  realtionship  any  S will 
yield  an  appropriate  when  solving  for  the  a^  in 
Equation  (32).  With  this,  both  equations  of  the  missile 
kinematics  are  completely  defined. 

Tracker  Noises . The  tracker  being  considered  in 
this  study  is  modeled  as  a typical  radar  for  a fighter  type 
aircraft.  A modeling  of  the  uncertainty  of  the  radar 
measurements  follows  a similar  development  as  was  described 


i 


for  the  glint  scintillation  and  thermal  noises  of  the 
missile  seeker.  The  strengths  and  correlation  times  for  the 
tracker  noises  are  given  in  Table  III. 

Table  III 

Sigmas  and  Time  Constants  for  Measurement  Noise 


Measurement 

Sigmas  ( a) 

Time  Constant 

6 (radians) 

uncorrelated  (Thermal) 
correlated  (Glint) 

0.00126 

0.00168 

0.5(sec) 

R (feet) 

uncorrelated  (Thermal) 
correlated  (Glint) 

11.7 

10.0 

0.5  (sec) 

• 

R (feet/sec) 

uncorrelated  (Thermal) 
correlated  (Glint) 

7.0 

4.242 

1.0  (sec) 

Guidance  Strip 

As  will  be  shown  in  Chapter  III,  the  major  thrust  will 
be  directed  at  modeling  the  "guidance  strip."  The  guidance 
strip  includes  the  seeker,  the  guidance  system  and  the  auto 
pilot/airframe  dynamics  of  the  missile.  It  is  the  set  of 
components  which  process  the  LOS  rate,  Q',  into  a developed 
lateral  acceleration.  Figure  13  shows  the  complete  truth 
model  guidance  strip. 

Summary 

This  chapter  presented  the  major  equations  and  assump- 
tions contained  in  the  truth  model  used  in  the  simulation. 
The  explicit  development  of  the  equations  has  been  omitted 


since  a de'tailed  derivation  can  be  found  in  Lutter  (Ref 
5:7-50) . Significant  detail  has  been  incorporated  into  the 
truth  model  to  yield  appropriate  results. 


i 
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III.  FILTEP>  KODELIPG 


Introd  ction 

The  !;asic  filter  model  considered  v.’itr;,Ln  this  study  v;os 
originally  developed  by  Lutter  in  Reference  5.  This  filter 
modeled  the  following  four  dynamic  states,  ;ind  t\.'0  parameters 
of  cin  air-to-air  missile. 

0,p  - LOS  angle  as  seen  by  the  ownship 

R - range 

R - range  rate  (relative  closing  velocity,  where 

by  convention:  R -V^  ) 

- x-velocity  of  missile  in  inertial  frame 
n - propoj'tional  navigation  constavt 

M/S  ■ “ ratio  of  the  mass  of  the  missile  to  crc's 
sectional  reference  area 

These  modeled  stator,  and  parameters  are  collectively  referred 
to  within  this  report  as  the  "basic  filoer,"  since  thf-y  serve' 
as  the  basis  for  all  other  designs  presented. 

A medeling  assumption  vrlthin  this  filter  ignortnl  the 
prefiitering  of  the  guidance  system  aJ.ong  v;ith  the  log  ef focus 
of  the  seeker  and  autopilot.  This  assumption  was  used  to 
dsvo-lop  the  simplest  filter  model  possible.  It  v;a3  argued, 
tl'.at  if  the  results  of  this  filter  were  unsatisfactory,  a 
more  refined  model  would  be  proposed.  The  assumption  essen- 
tially equated  lateral  acceleration  to  commanded  acceleration. 


Recalling  the  proportional  navigation  law,  this  assumption 


is  e:<pros3cd  in  the  following  equation. 


IBIS  FACE  IS 
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= nV  0 
c 


(36) 


As  previously  discussed  in  Chapter  2,  the  term  "guidance 
strip"  refers  to  those  components  which  process  LOS  information, 
(6'),  into  developed  lateral  acceleration.  Figure  14  shows 
the  actual  guidance  strip  in  the  missile  along  with  the  model 
based  upon  Equation  (36)  and  used  by  Lutter  in  Reference  5. 

The  effects  of  aberration  error  and  seeker  noise  upon  6 in 
the  truth  model  are  indicated  by  the  use  a prime.  Also,  the 
subscript  "s"  has  been  used  to  indicate  the  effects  of  the 
seeker  dynamics  upon  9 ' . 

The  analysis  within  this  report  continues  to  use  the 
basic  filter  developed  by  Lutter  with  various  models  of  the 
guidance  strip  inserted  into  the  filter.  For  purposes  of 
identification,  these  various  models  of  the  guidance  strip 
are  used  to  define  the  "order  of  the  missile,"  which  is 
modeled  in  the  filter.  The  number  of  poles  in  the  guidance 
strip  transfer  function  determines  the  missile  order.  In 
this  way,  the  basic  filter  designed  by  Lutter  is  called  the 
"zero  order  missile  filter." 

The  performance  of  the  zero  order  filter  was  considered 

inadequate  for  accurate  pointing  and  tracking,  which  was  the 

original  design  goal  for  the  filter.  Figure  15  shows  the 

estimation  of  the  state  for  the  zero  order  filter  when 

mx 

it  was  given  perfect  knowledge  of  n and  M/S.  In  general, 
knowledge  of  these  parameters  will  not  be  available  and  the 
filter  will  have  to  estimate  their  value.  Thus,  Figure  15 
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Zero  Order  Missile  Filter  Guidance  Strip 


I 


is  portraying  the  best  estimation  possible  for  the  zero  order 
filter  when  the  missile  flies  a high-g  trajectory.  The  box 
in  the  upper  right  hand  corner  of  the  figure  contains  the 
important  information  about  the  curve.  For  this  plot,  it  indi- 
cates that  the  center  curve  is  the  "mean  error"  of  the 

nuc 

state  (x- velocity)  when  all  parameters  were  set  equal  to 
their  true  values.  It  also  indicates  the  use  of  data  set 
120  for  the  zero  order  missile  filter.  The  data  sets  are 
explained  in  Appendix  B.  The  center  curve  is  an  ensemble 
average  of  the  mean  error  made  by  the  filter  over  20  simu- 
lations of  the  high-g  intercept  scenario.  The  mean  error 
is  enveloped  by  plus  and  minus  one  standard  deviation.  Thus, 
fora  time  of  approximately  1.7  seconds  in  the  time  history 

of  this  curve,  the  filter  error  in  the  state  can  be 

mx 

observed  to  be  150  ft/sec  with  a one  sigma  confidence  of 

25  ft/sec.  The  x-velocity  state,  , of  each  filter 

mx 

designed  in  this  report  was  the  most  sensitive  state  to 
filter  changes.  For  this  reason,  it  will  be  used  to 
represent  the  overall  performance  of  the  various  filters 
to  be  discussed. 

It  was  suggested  by  Lutter  that  modeling  the  lag 
effects  of  the  autopilot  would  improve  the  estimates 
achievable  by  the  zero  order  filter.  This  conclusion  was 
based  upon  the  observation  that  the  filter  was  estimating 
too  large  an  initial  deceleration  value  of  the  missile  in 
the  - direction,  causing  an  initial  positive  bias.  This 
condition  occurred  because  of  the  zero  delay  for  lateral 


acceleration  modeled  in  the  filter  (Ref  5:84-85) . 


The  initial  modeling  effort  in  this  report  was  aimed  at 
acquiring  a proper  benchmark.  The  benchmark  would  include 
as  complex  a model  as  necessary  to  provide  accurate  state 
estimation.  In  this  way  the  benchmark  would  represent  a 
lower  bound  for  filter  error  and  would  be  used  as  a com- 
parison when  evaluating  reduced  order  filters.  The  approach 
taken  in  designing  the  benchmark  was  based  upon  the  sugges- 
tion made  by  Lutter,  but  it  encoximpassed  a more  exhaustive 
set  of  possible  models  for  the  guidance  strip.  Figure  16 
represents  the  possible  models  th<  ; were  considered  for 
each  component  of  the  guidance  strip.  A column  labeled 
"a^  transferred"  is  also  included  in  this  figure.  This 
column  was  used  only  in  conjunction  with  an  autopilot  model 
and  refers  to  an  ad  hoc\  procedure  in  which  the  truth  model 

value  of  commanded  acceleration,  a , is  passed  directly  to 

\ ^ 

the  autopilot  of  the  filter.  This  procedure  was  used  in 
designing  a benchmark  and\  will  be  discussed  in  depth  later 
in  this  chapter.  Combining  a model  choice  from  each  compo- 
nent column  or  choosing  an\  autopilot  model  with  commanded 
acceleration  transferred,  defines  a proper  benchmark  design. 

Once  a benchmark  had  been  decided  upon,  the  insights 
gained  from  its  performance  inspired  the  development  of  two 
filters,  a fourth  order  missile  filter  and  a first  order 
missile  filter.  The  fourth  order  missile  filter  included 
exact  models  of  the  guidance  system  and  autopilot.  This 
filter  was  used  to  determine  filter  sensitivity  to  the 
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navigation  constant,  n,  the  M/S  ratio,  and  the  bandwidth  of 
the  guidance  strip.  It  also  established  the  best  estimation 
possible  when  given  perfect  knowledge  of  these  parameters. 

The  first  order  filter  was  proposed  as  a reduced  order 
filter  which  would  provide  accurate  estimation  of  the 
missile's  states  and  three  of  its  key  parameters.  This 
filter  used  a first  order  lag  to  approximate  the  complete 
guidance  strip  of  the  missile.  The  performance  of  the  first 
order  filter  was  unexpectedly  impressive,  considering  that 
the  first  order  missile  model  in  the  filter  was  approximating 
the  entire  fifth  order  guidance  strip.  A complete  discussion 
as  to  why  this  filter  performed  so  well  will  be  deferred  un- 
til Chapter  IV;  however,  a foundation  for  its  understanding 
can  be  gained  from  the  modeling  presentation  given  in  this 
chapter . 

Figure  17  presents  the  guidance  strips  of  all  filters 
that  were  analyzed  in  this  report.  The  guidance  strips  are 
arranged  in  order  of  complexity  starting  with  the  zero  order 
filter.  This  chapter  will  discuss  the  rationale  for  each 
design  in  chronological  order  of  their  development,  after 
reviewing  the  modeling  of  the  basic  filter. 

Modeling  of  the  Basic  Filter 

The  use  of  the  Extended  Kalman  Filter  algorithm  requires  | 

that  the  system  of  interest  be  described  by  a set  of  first 
order  ordinary  differential  equations  and  that  the  random 


disturbances,  which  model  the  system  uncertainty,  be  input 
in  a linear  additive  fashion  (Ref  7:179).  Equation  (37) 


Truth 

Model 


V 

n c 


a-s  +a,s+a 


! 

I 


i 
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expresses  this  modeling  approach  in  vector  notation 

x(t)  = f (x(t)  ,u(t)  ,t)  + G(t)w(t)  (37) 

where 

x(t)  = n-dimensional  state  vector  (the  underscoring 
indicates  that  the  variable  is  a vector) 
w(t)  = zero  mean  white  Gaussian  noise 

with 

E {w(t)w'^(t  + T)}  = Q(t)5(T) 

Restricting  the  modeling  for  the  filter  to  this  general  form, 
the  following  subparagraphs  develop  and  define  the  first 
six  states. 

0 

Line  of  Sight  Angle  ( T) . Figure  18  presents  the 
inertial  frame  with  a depiction  of  the  ownship  and  missile 
at  some  time  after  the  start  of  the  intercept.  In  this 
figure,  the  position  of  the  ownship  and  missile  are  defined 
by  the  two  inertial  vectors  Rq,j,,  and  respectively.  The 
vector  R',  is  aligned  with  the  x-axis  of  the  ownship  antenna 
and  is  the  range  of  the  missile  from  the  ownship.  Using 
vector  addition 


-OM  -OT  — 


(38) 


Utilizing  the  theorem  of  Coriolis,  and  notation  established 
in  reference  14,  the  inertial  velocity  of  the  missile  is 
found  to  be 


“ ^TLOS^OM  “l(TLOS)  * ^qm 


(39) 
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After  siubstituting  Eq’iation  (38)  into  (39)  and  reducing, 
Eg'-ation  (39)  can  be  written  as 


PR  =PR  H-P  R + u,  xR 

lk)M  I^T  ‘TLOG-  -I(TLOS)  -- 


or  equivalently , coordi;uiti;:ed  in  the  I frarae, 


+ s'  + <Si(Tic3)  s) 


where 


P 

TLOS 


_ time  derivative  with  respect  to 
the  inertial  frarri-e 

= tiiTiG  derivative  with  respect  to  the 


TLOS  frame 


= relative  closing  velocity  between 
missile  and  ownship 


- inertial  velocity  vector  of  the  missile 


= inertial  velocity  vector  of  the  tracker 


-I  (TLOS; 


= tile  angv.  lar  rate  bet'ween  the  I frame 


and  the  TLOS  frame 

Expressing  Lauation  (-11)  in  terms  of  its  components, 
and  solving  for  the  magnitude  of  j f follciwing 
two  equal  ions  car  be  found. 


“l(TLOS)  = § cote  + {42a) 

K P K 

O'  (TLOS)  ---J^tanl,,  H-  (42b) 

i i T> 
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where  by  definition  (see  Figure  18) 


the  deoirod  modeling  of  LOS  angle  has  been  achieved.  It  was 
decided  to  use  Equation  (42a)  when  sin  ^ .707  and  Equa- 
tion (42b)  v.'hen  sin  0^  <.707  to  avoid  the  singulari ties  of 


Equation.^  (42a)  and  (42b)  are  also  used  to  define  an 
important  relationship  which  will  be  used  to  solve  for  the 


y-componen t of  the  missil.e  velocity  in  the  inertial  front!. 
By  setting  Equ<. Lions  (42a)  and  (42b)  equal  to  one  another, 

, J 

solving  for  V ' yields 
my 


R CObd, 


CO30_,  -V  cosS^  + RSine, 
mx T_ 

SinG„ 


X 

Equativon  (44)  shows  that  V car,  ba  found  dircctlv  v.'ith 
^ my 

I 

linowlcdgo  of  V , 6„,  and  R,  This  ecuation  is  used  outside 
,mx ' T 

the  structure  of  the  filter  and  thus  eliminates  the  need  for 
estimating  both  components  of  the  inertial  velocity  of  the 
missile.  The  modeiing  in  this  analysis  assumed  that  no 
ccndit''.ons  exist  v.-hioh  W'Ould  favor  estimating  one  component 
of  inerti.al  velocity  over  the  otner.  The  x-coiviponent  of 
the  velocity  was  cscinated  in  all  cases. 

X-component  ^ 

ing  as  a first  order  d.-f  f crenti  a i equation  c,an  be 
accomplished  by  resolving  the  acceleration  acting  upon  the 
missile  into  its  x and  y cor.'ponents , figure  19  shows  ft 


1 

I 


( 
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Diagram  gj:  the  MissiB.e 


free  Lo  y diagran  of  the  nissile  drawn  in  the  inertial  frar.e. 

• I • T 

Solving  for  V and  V 

m:<  my , 


V - -a^  Sin  Y - a,  cosy 
h’X  D ' ni  L ' in 


(45a) 

{45b) 


where 

aj^  = total  drag  acting  upon  the  n.issile 

= lateral  acceleration  due  to  lift 
Ju 

Y = angle  defining  incidence  of  inertial  velocity  of 
missile  to  zc-axis  of  inertial  frame. 

It  is  important  to  erap' asize  here  that  the  modeling  assump- 
tions introduced  in  Chapter  I allow  for  specification  of 
the  acceleration  terms  in  Equations  (45a)  and  (45b). 
Repeated  again,  these  assumptions  are: 

1)  that  the  missile  v;ill  be  using  proportional  naviga- 
tion 
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2)  that  the  missile  is  non- thrusting. 

* 1 * I 

Equation  for  V is  developed  the  same  as  V for  use  in 
^ my  mx 

defining  the  range  rate  state  in  the  following  section. 

• T 

Proceeding  with  the  development  of  an  expression  for 

a^j  is  determined  starting  with  the  basic  aerodynamic 
equation  for  total  drag. 

The  use  of  Equation  (46)  is  dependent  upon  knowledge  of 
(l/CM/S))  and  C^.  The  other  terms  can  be  considered  known  or 
can  be  obtained  from  measurements  of  the  missile's  kinematics 
The  term  (M/S)  is  given  special  attention  in  a later  discus- 
sion and  for  the  moment  will  be  considered  known.  Thus 
the  coefficient  of  drag,  Cj^,  remains  to  be  defined  in 
terms  which  are  observable  to  the  ownship  either  directly 
or  through  some  filtering  process.  This  approach  is 
assuming  that  apriori  knowledge  of  Cj^  is  not  available. 

The  coefficient  of  drag  can  be  expressed  as  the  sum  of 
the  zero  lift  drag  and  the  induced  drag  coefficients. 

S " ^Do  Si  (47) 

where 

CjjQ  = coefficient  of  zero-lift  drag 

Cjjj  = coefficient  of  induced  drag 

The  approximation  chosen  to  represent  C_^  and  C are 

UU  ul 
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where 


A = speed  of  sound  at  the  appropriate  altitude  (Ref  5:63) 


2 2 
p S A 


Substituting  in  the  total  expression  for  Cj^  into  Equation 
(46)  and  rearranging  terms  yields 


a^  = JspS/M  r2A’*  , 4M^ 

° “212 

P S 


V^2 

p^S^A 

m 


In  the  above  equation  p and  A are  considered  known  and  for 
this  analysis  are  held  constant  for  an  altitude  of  10,000 
feet  above  sea  level  for  a clear,  standard  day.  Incorporating 
the  complete  expression  for  a^^  into  Equations  (45a)  and  (45b) 
yields  the  desired  relationship. 


‘ <!'  * 'Ia*  'f>  -®L  <51) 

m 


The  modeling  for  a^  in  Equation  (51)  is  described  separately 

for  each  filter  since  the  definition  of  this  state  will  vary 

cimong  the  different  models  of  the  guidance  strip. 

The  total  inertial  velocity  of  the  missile  and  the 

angle  Y are  calculated  with  the  current  estimates  of 
^ xci  inx 

and  by  the  following  equations: 

+ (52] 

m \ mx  my 
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tan-^  (^) 

V 


The  term  M/S  is  an  unknown  parameter  in  Equation  (51) 
and  is  estimated  by  the  filter.  The  modeling  for  this  term 
will  be  discussed  later  in  this  chapter. 

Modeling  for  Range  and  Range  Rate  (R, R) . The  modeling 
for  range  rate  is  obtained  by  taking  the  derivative  of 
Equation  (40) , (repeated  here) 


PAm  ^I^T  ^TLOS-  -I(TLOS)^  - 


with  respect  to  the  I-frame.  Applying  the  Theorem  of  Coriolis 
Equation  (54)  is  derived 


Pi^Som  “ Pi  5ot  PtLOS-  ^-I(TLOS)  * PtLOS- 


PtLOS  -KTLOS)  “l(TLOS)  * (TLOS)  ^ 

Equation  (54)  can  be  cc ^rdinatized  in  the  I-frame  and  re- 
%n:itten  as 

5^  * B ■*■2  ^-i(TLos)  ^Ptlos  -I(TLOS)  * 


* ^-I(TLOS)  ^ ^-I(TLOS)  B))  ( 


where 


= inerbial  missile  acceleration 
—in 


a^  = inertial  tracker  acceleration 

•*T 

R = relative  acceleration  along  the  LOS 
-I(TLOS)  * 5 ” Coriolis  acceleration  term 

^TLOS  -I(TLOS)  * B = tangential  acceleration  term 
-I(TLOS)  * ^-I(TLOS)  ^ ^ centripetal  acceleration  term 

The  Coriolis  and  tangential  acceleration  components  are 

TTjOS 

perpendicular  to  the  X direction  and  do  not  contribute 

TIjOS 

to  the  X acceleration  component.  Equation  (55)  then 
reduces  to: 


f + (“ktlos)  * ‘Sktlos)  * S'i'  (5 

• • 

Coordinatizing  in  the  TLOS  frame,  and  solving  for  R (accel- 

TTjOS 

eration  along  the  X axis) , the  desired  equation  for 
modeling  range  rate  is  derived. 


R = 


2 

<0  R 


+ a 


TLOS 

mx 


^TLOS 

^Tx 


(57) 


where 

“ “ -I (TLOS) 

TLOS 

The  acceleration  term  a^^  is  considered  known  through  a 
coordinate  transformation  of  the  accelerations  obtained  from 


the  inertial  navigation  unit  on  board  the  ownship.  The 
accelerations  obtained  from  the  navigation  unit  are  con- 
sidered deterministic  since  the  precision  of  these  measure- 
ments is  much  greater  than  that  available  from  a typical 

TLOS 

aircraft  radar,  a is  found  through  the  use  of  a direction 

mx  ^ 

* I *1 

cosine  matrix  transforming  V and  V into  the  TLOS  frame. 

mx  my 

Recalling  that  9^  is  a measure  of  the  angle  between  the  TLOS 
and  inertial  frame,  the  proper  transformation  is 


TLOS  ^ cos0„  + sine  (58) 

mx  mx  T my  T 

* I * I 

Using  the  expressions  for  V_  and  V found  in  the  previous 

mx  my 

TLOS 

section,  and  substituting  the  complete  expression  of  a^^ 
back  into  Equation  (57) , the  desired  state  modeling  of  range 
rate  is  obtained; 

3 ^2 

R = e^R  + [{pA**  (S/M)  siny^  - (|^)  (M/S)  ^ siny^^^ 


1 

- a^^  cosy^jj}  cos0^  + {-  pA^  (S/M) 


(— ) (M/S)  cosy^  - a^^  siny^}  sine,^] 


TLOS 


Equation  (59)  is  modeled  in  the  filter  as  two  first  degree 

equations : 

where 


5t 


and  Proportional  Navigation  Constant  (M/S,  n) . This  study 
has  considered  both  the  M/S  ratio  and  the  navigation  con- 
stant to  remain  invariant  throughout  the  scenario.  For  this 
reason  these  two  parameters  could  be  modeled  as  random 
biases. 

(M/S)  (t)  = 0 (62] 

n(t)  = 0 (63] 

Once  the  filter  has  obtained  an  initial  estimate  of  the 
constant,  the  random  bias  model  causes  the  Kalman  gain  to 
approach  zero  and  essentially  ignore  available  information 
brought  in  after  this  time.  Since  the  complete  missile 
model  in  the  filter  is  not  exact,  it  is  possible  for  the 
filter  to  acquire  an  incorrect  initial  estimate  of  these 
constants.  This  behavior,  referred  to  as  "learning  the 
wrong  state  too  well,"  can  be  avoided  by  adding  a pseudo 
noise  to  each  channel  and  thus  avoid  the  condition  of  having 
the  Kalman  gain  go  to  zero  on  these  channels  (Ref  3:189; 
6:204).  Thus,  the  complete  model  of  these  parameters  is 


defined  as 


(64) 


(M/S) (t)  = w^(t) 
n (t)  = W2(t)  (65) 

where 

E [Wj^(t)  (t+T)]  = Qj^5(t) 

E [W2(t)  W2’^  (t+T)]  = Q26(t) 

with  and  suitably  chosen  during  the  tuning  process. 

Searching  for  a Benchmark.  Figure  16,  presented  in  the 
introduction  of  this  Chapter,  laid  out  the  possible  alter- 
natives for  designing  a benchmark.  The  intent  in  finding  a 
benchmark  was  to  demonstrate  the  best  performance  obtainable 
from  even  the  most  complicated  form  of  a potential  filter 
design.  The  benchmark  is  sought  not  for  possible  implementa- 
tion, but  more  for  a comparison  between  it  and  reduced  order 
filters.  In  some  cases,  as  in  this  analysis,  the  benchmark  is 
given  knowledge  of  a system's  structure  along  with  its  exact 
parameter  values,  all  of  which  would  not  be  available  should 
its  design  be  considered  for  implementation.  This  image  of 
a benchmark  is  what  prompted  the  development  of  the  three 
designs  that  follow.  The  final  design  was  chosen  to  represent 
the  benchmark  for  all  possible  filter  designs  that  could  be 
implemented . 
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Fig.  20.  Second  Order  Missile  Filter 


Second  Order  Design.  The  first  design  considered  lag 
characteristics  of  the  autopilot  as  the  most  important  effect 
influencing  filter  performance.  Although  the  actual  autopilot 
was  third  order,  it  was  adaptively  controlled  by  feedback 
gains  (refer  to  Chapter  II)  to  respond  as  a second  order 
system  with  an  effective  damping  ratio  of  .707  and  a natural 
frequency  of  7.07  rad/sec.  These  values  are  based  upon 
previous  design  work  for  this  missile  model  (Ref  5:26).  The 
choice  of  the  second  order  model  was,  therefore,  a reasonable 
first  approach  in  capturing  the  lag  effects  of  the  autopilot 
(see  Figure  20)  . 

The  performance,  however,  was  marginal.  The  estimates  in 
the  state  were  greater  then  100  ft/sec.  This  could 

hardly  be  used  as  a benchmark  since  this  would  be  admitting 
that  any  filter  designed  for  implementation  would  in  the  best 
of  circumstances  do  no  better  then  100  ft/sec. 
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Fig.  21.  Third  Order  Missile  Filter 

Third  Order  Design.  The  second  attempt  in  defining  a 
benchmark,  used  the  exact  structure  of  the  autopilot.  This 
approach  still  considered  the  lag  effects  of  the  autopilot 
as  the  essential  effect  to  be  modeled.  The  difference  between 
this  approach  and  the  second  order  design  was  that  the  actual 
autopilot  was  nonminimum  phase.  This  is  usually  true  of  tail- 
controlled  missiles  since  the  initial  movement  of  the  rear 
fins  causes  acceleration  in  the  opposite  direction  desired. 
This  "tail-wags-the-dog"  characteristic  shows  up  in  the  auto- 
pilot transfer  function  as  a right-half-plane  zero.  It 
seemed  possible  that  inclusion  of  this  effect  was  needed  for 
producing  a benchmark  (Ref  11:20). 

The  coefficients  of  the  autopilot  transfer  function  were 
time  varying  since  it  was  designed  to  be  adaptive.  An 
ensemble  average  of  each  coefficient  for  six  tine  points  of 
the  high-g  intercept  was  made  from  twenty  runs  of  the 
simulation.  Table  IV  shows  the  poles  and  zeroes  of  the 
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Table  IV 
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sec  -47.77  -16.69  -4.06  -35.99  33.33 


transfer  function  for  the  six  time  points.  Even  though  the 
coefficients  did  vary,  the  actual  dominant  pole  and  zero 
movement  was  slight  enough  that  one  set  of  coefficients  could 
well  approximate  the  autopilot  for  the  entire  scenario. 
Although  this  was  true  for  steady  state  estimates,  it  was 
found  that  using  the  coefficients  for  t = 0 slightly  reduced 
the  transients  in  the  state  error  plots.  These  coefficients 
were  therefore  used.  The  nonminimum  phase  characteristic  of 
the  autopilot  was  modeled  by  using  constant  coefficients  and, 
the  ease  of  implementing  this  approach  made  it  more  attractive 
compared  to  the  integration  of  the  complete  adaptive  auto- 
pilot (see  Figure  21) . 

Quite  unexpectedly,  the  third  order  design  performance 

was  not  much  improved  over  the  second  order  design.  The 

mean  error  in  channel  was  still  reaching  a 110  ft/sec 

maximum  (all  plots  for  this  design  can  be  found  in  Appendix 

B) . Again  this  design  was  considered  unacceptable  for  a 

benchmark.  Although  no  minimum  error  had  previously  been 

decided  upon,  the  fundamental  structure  of  the  filter  was 

considered  sound  and  capable  of  producing  better  estimates 

of  the  state, 
mx 

Benchmark.  Proceeding  with  the  assumption  that  the  other 
modeling  in  the  basic  filter  was  adequate,  the  true  a^  of  the 
missile  was  passed  to  the  filter  design  which  modeled  the 
third  order  autopilot.  This  was  done  in  an  attempt  to 
determine  if  more  modeling  of  the  guidance  strip,  prior  to 


the  autopilot,  was  necessary.  Also,  if  this  approach  failed 
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Fig.  22.  Benchmark 


to  produced  reasonably  accurate  estimates,  it  would  indicate 

a major  flaw  in  the  fundcimental  design  structure  of  the  filter. 

The  benchmark  is  shown  in  Figure  22. 

The  performance  improvement  of  this  design  was  dramatic 

when  compared  to  the  second  and  third  order  designs.  For 

exaunple,  the  mean  error  committed  in  the  channel,  the 

mx 

most  sensitive  state,  did  not  exceed  10  ft/sec.  This 
enhancement  in  performance  confirmed  the  idea  that  more 
modeling  of  the  guidance  strip  was  necessary.  It  also 
indicated  that  the  autopilot  was  not  the  critical  subsystem 
of  the  guidance  strip  and  that  a simpler  autopilot  model 
would  suffice. 

Although  this  ad  hoc  design  could  never  be  implemented, 
it  was  considered  to  be  a proper  benchmark  and  would  be  used 
as  a comparison  for  reduced  order  filters.  The  small  errors 
committed  in  each  channel  would  serve  as  a critical  test  on 


Fig.  23.  Fourth  Order  Missile  Filter 


the  performance  of  all  filters  analyzed  for  implementation. 

The  complete  set  of  benchmark  plots  are  contained  in  Appendix 

B. 

Fourth  Order  Missile  Model . The  benchmark  design  provided 
a lower  bound  for  filter  errors  and  indicated  the  need  for 
more  modeling  of  the  guidance  strip  prior  to  the  autopilot. 

The  ad  hoc  approach,  however,  used  in  this  development  did 
not  include  models  of  the  seeker  and  guidance  system.  The 
development  of  implementable  reduced  order  filters  would  need 
specific  knowledge  of  the  dominant  effects  of  these  two 
components.  To  gain  this  insight,  the  fourth  order  missile 
model  was  proposed. 

Upon  examination  of  the  three  components  of  the  guidance 
strip,  it  was  observed  that  the  guidance  system  had  a limiting 
bandwidth  of  3 rad/sec.  Comparing  this  to  the  7 rad/sec  band- 
width of  the  autopilot  and  the  50  rad/sec  bandwidth  of  the 
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Pig.  24.  First  Order  Missile  Filter 

seeker,  it  was  hypothesized  that  the  modeling  of  the  guidance 
strip  should  include  the  bandlimiting  characteristic  of  the 
guidance  system.  The  guidance  strip  of  this  filter  model  is 
shown  in  Figure  23. 

The  fourth  order  missile  filter  produced  estimates  of  the 
state  with  a maximum  mean  error  of  25  ft/sec.  This  was 
only  15. ft/sec  greater  then  the  benchmark  filter,  in  which 
a^  was  considered  to  be  known  exactly.  This  performance 
emphasized  the  need  to  model  the  bandlimit  effect  of  the 
guidance  system. 

First  Order  Missile  Filter.  The  significant  performance 
of  the  fourth  order  filter  suggested  the  possibility  of  a 
reduced  order  filter  which  modeled  the  complete  guidance 
strip  as  a first  order  lag.  This  reasoning  followed  since 
the  guidance  system,  itself  a first  order  lag,  appeared  to  be 
the  dominant  component  of  the  entire  guidance  strip.  It  was, 
therefore,  decided  to  approximate  the  complete  guidance  strip 


as  a first  order  lag  (Figure  24),  With  this  approach,  a fifth 
order  guidance  strip,  considering  time  invariant  coefficients, 
was  modeled  by  a first  order  lag  in  the  filter. 

An  iterative  search  was  performed  to  determine  a time 
constant,  t^,  which  would  produce  the  least  error  in  the 
state  estimates.  During  this  search,  the  proportional 
navigation  constant,  n,  and  the  M/S  ratio  were  assumed 
perfectly  known.  The  results  of  this  search  showed  that 
the  best  estimates  were  obtained  with  a time  constant  of 
.85  seconds.  The  maximum  error  of  the  first  order  missile 
model  was  35  ft/sec,  only  10  ft/sec  greater  than  the  fourth 
order  filter  and  25  ft/sec  greater  than  the  benchmark. 

This  filter  was  the  first  attempt  at  providing  a reduced 
order  model  of  the  guidance  strip  which  could  possibly  be 
implemented.  However,  the  design  required  knowledge  of 
for  the  particular  missile  being  tracked.  It  also  needed 
values  for  n and  M/S.  If  the  filter  was  insensitive  to  these 
parameters,  it  might  be  possible  to  use  nominal  values  for 
all  three.  Of  course,  if  the  filter  was  highly  sensitive  to 
some  or  all  of  these  paraimeters,  their  estimation  might 
improve  the  overall  performance  of  the  filter.  For  this 
case,  a random  walk  would  be  chosen  to  model  Tj.  This  model 
would  portray  the  essential  characteristic  of  t^,  which  is 
assumed  constant  over  the  entire  scenario.  As  with  n and 
M/S  the  pseudo  noise  is  added  to  avoid  having  the  Kalman  gain 
components  affecting  the  estimates  go  to  zero.  The 

complete  model  is  shown  below. 


Tf  (t)  = W3(t)  (66)  ^ 

where  E [w^Ct)  W3(t+T)]  = Q^fiCr) 

The  value  of  would  be  chosen  during  tuning  of  this  filter.  j 

Summary.  This  Chapter  presented  the  derivation  of  the  i 

i 

equations  for  the  first  six  states  common  to  all  filters.  A 

benchmark  filter  was  designed  which  gave  insight  into  the  j 

important  areas  of  modeling.  The  development  of  a fourth  i 

4 

order  missile  filter  resulted  from  this  insight.  In  an 
attempt  to  reduce  the  complexity  of  the  filter's  guidance 

j 

strip  model,  a first  order  lag  was  chosen  to  approximate  the  I 

i 

complete  guidance  strip  of  the  missile.  j 

I 

The  time  constant,  t^,  of  the  first  order  lag  was  chosen 
to  minimize  the  error  in  the  state  estimates.  The  design, 
however,  required  a detailed  analysis  concerning  the  three 
parameters,  n,  M/S,  and  to  completely  evaluate  its 
performance.  A sensitivity  analysis  for  these  uncertain 

i 

parauneters  could  possibly  lead  to  an  adaptive  filter  that 
estimates  some  (most  important)  or  all  (depending  upon  the 
computational  resources  available)  of  the  uncertain 

parameters.  j 
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IV.  Sensitivity  and  Parameter  Estimation  Analysis 


Introduction 

The  modeling  effort  discussed  in  Chapter  III  led  to  two 
filters  of  practical  interest;  the  fourth  order  missile 
filter  and  the  first  order  missile  filter.  The  fourth  order 
filter  demonstrated  the  benefits  of  including  a model  of  the 
guidance  system  and  a detailed  model  of  the  missile's  auto- 
pilot. This  extensive  modeling  required  perfect  knowledge 
of  many  parameters  which  in  most  cases  are  uncertain. 
Implementation  of  such  a filter  is,  therefore,  considered 
impractical.  However,  this  detailed  modeling  would  be 
advantageous  in  determining  the  sensitivity  of  the  filter 
design  to  variations  in  key  parameters  of  the  missile. 

The  first  order  filter  is  more  attractive  for  implemen- 
tation since  it  approximates  the  complete  guidance  strip  as  a 
first  order  lag.  It  does  not  require  specific  knowledge  of 
the  structure  and  parameters  of  the  particular  autopilot,  or 
the  time  constant  of  the  guidance  system.  Instead,  the 
significant  aspects  of  the  missile  are  modeled  with  only 
three  uncertain  pareimeters: 

1)  n - the  proportional  navigation  constant 

2)  M/S  - the  ratio  of  mass  to  cross  sectional 

area 

3)  - the  time  constant  of  the  first  order  lag. 
It  should  be  remembered  that  a fourth  uncertain  parameter, 

Cjj  was  redefined  in  terms  of  quantities  that  could  be  measured 
or  estimated  (refer  to  Chapter  III) . 
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The  design  approach  in  this  analysis  considered  two 
possible  alternatives  for  determining  each  unknown  parameter. 
Either  a nominal  value  would  be  chosen  for  the  parameter, 
or  the  filter  would  perform  on-line  estimation  to  acquire  the 
value  of  the  parameter.  The  criterion  used  to  decide  this 
issue  was  the  set  of  results  of  a sensitivity  analysis  per- 
formed on  the  fourth  order  filter. 

This  study  concluded,  from  the  results  of  this  sensitiv- 
ity analysis,  that  on-line  estimation  of  each  parameter  was 
both  possible  and  beneficial.  A systematic  approach  was 
then  devised  which  gradually  increased  the  scope  of  the 
estimation  problem.  This  approach  consisted  of  a sequence 
of  estimation  problems  made  more  complex  by  the  inclusion  of 
additional  parameters  to  be  estimated.  Eventually,  this  led 
to  the  simultaneous  estimation  of  all  five  states  and  all 
three  parameters  of  the  first  order  missile. 

Chapter  IV  begins  with  a discussion  on  the  tuning 
philosophy  used  throughout  this  study.  This  is  followed  by 
a presentation  which  justifies  the  use  of  twenty  simulation 
runs  for  the  Monte  Carlo  analysis.  The  Monte  Carlo  analysis 
was  used  to  determine  the  performance  of  each  filter  design 
analyzed.  Following  this  is  a detailed  explanation  of  the 
sensitivity  analysis  of  the  fourth  order  filter.  The  final 
topic  covers  the  parameter  estimation  analysis  of  the 
first  order  filter. 

Tuning  Philosophy 

The  tuning  philosophy  used  throughout  this  study  can 
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be  expressed  as  "covariance  matching".  Basically,  this  is 
the  matching  of  the  covariance  history  generated  by  the  filter 
from  the  Kalman  algorithm  to  the  variances  in  the  filter 
errors  which  are  calculated  in  the  Monte  Carlo  analysis. 

The  covariance  in  the  Kalman  algorithm  is  an  indication  of 
the  confidence  of  the  filter's  estimate.  Covariance  match- 
ing, then,  tunes  the  covariance  of  the  filter  to  reflect  the 
true  confidence  in  its  estimates. 

The  "tuning"  of  the  filter  was  performed  by  a proper 
choice  of  pseudo  noise  strength,  Q,  and  the  initial  filter 
covariance,  Pq,  for  each  state.  The  initial  value  of 
measurement  noise  strength,  modeled  by  R,  was  derived  from  a 
priori  knowledge  of  the  ownship's  radar  performance.  R was 
also  considred  a tuning  parameter,  but  its  initial  value, 
described  above,  was  used  as  a lower  bound.  This  was 
necessary  to  avoid  the  condition  in  which  the  filter  weighted 
the  incoming  measurements  with  more  certainty  than  the  radar, 
as  modeled,  was  capable  of  providing.  These  tuning  param- 
eters, as  well  as  the  initial  state  estimates,  are  included 
with  their  associated  error  plots  in  Appendix  B. 

During  the  tuning  of  the  filter  for  parameter  estimation, 
it  became  obvious  that  strict  adherence  to  the  philosophy  of 
covariance  matching  was  not  advantageous  for  the  parameter 
estimates.  The  pseudo  noise  strength,  Q,  needed  for 
adequate  parameter  estimation  was  too  large  for  matching  the 
filter  covariance  to  the  true  variance.  This  resulted  in  the 
filter  covariance  always  being  larger  than  the  true  variance 
for  the  parameter  estimates. 
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However,  this  insured  the  conservative  nature  of  the  filter  i 

by  never  attributing  more  confidence  to  its  estimates  than  the 
filter  was  actually  achieving,  thereby  reducing  the  risk  of 
divergence.  Since  the  tuning  was  performed  for  one  trajectory, 
adaptive  self  tuning  may  be  required  for  an  operational 
system  to  obtain  performance  comparable  to  the  results  pre- 
sented in  this  study. 

Variance  Convergence 

The  use  of  the  Monte  Carlo  analysis  requires  the  deter- 
mination of  an  appropriate  number  of  simulations  such  that 
the  statistics  presenting  filter  performance  are  generated 
with  a reasonable  confidence.  The  number  of  simulations  was 
initially  chosen  as  twenty.  It  was  decided  that  this  number 
would  be  increased  if  twenty  simulations  could  not  provide 
valid  results. 

To  confirm  that  twenty  simulations  were  adequate,  the 
variance  of  the  true  error  was  plotted  versus  the  number  of 
simulation  runs.  With  a proper  choice  of  simulation  runs, 
this  plot  would  show  little  change  in  the  variance  with  the 
addition  of  each  successive  run.  Figure  25  illustrates  such 
a plot  for  four  time  points  of  the  velocity  variance 
history.  This  plot  clearly  indicates  that  twenty  simula- 
tions were  an  adequate  choice  for  this  problem.  The 
addition  of  successive  runs  over  twenty  may  mathematically 
increase  the  confidence  of  the  statistics  but,  it  is 
obvious  that  there  will  be  little  change  in  the  value  of  the 
variance.  Similar  plots  of  other  states  can  be  found  in 


vmifwce 


4.0Q  7.0a  IQ.OQ  13.00  18.00  19.00  22.00 

RUN  


r 


Appendix  B which  also  illustrates  this  point. 

Sensitivity  Analysis 

The  fourth  order  filter  was  used  for  the  sensitivity 
analysis  since  the  detailed  modeling,  included  in  its  design 
of  the  guidance  strip,  allowed  for  an  easier  interpretation 
of  the  filter  performance.  All  states  were  initialized  with 
zero  error  and  the  high-g  trajectory  was  chosen  for  the 
simulation.  It  was  suspected,  and  later  observed,  that 
the  filter  would  be  less  sensitive  to  parameter  variations 
over  the  low-g  trajectory  (discussed  further  in  the  section 


on  estimation  of  n and  . The  three  parameters,  n,  M/S, 


and  1 2!  were  set  in  the  filter  to  the  original  values  used 
in  the  truth  model.  The  truth  model  values  were  then  sepa- 
rately varied  over  what  was  considered  to  be  their  maximum 
range.  The  fourth  order  filter's  sensitivity  to  variations 


in  ^2  would  also  be  an  indication  of  the  first  order  filter's 


sensitivity  to  t^,  while  n,  and  M/S  were  modeled  directly 
in  this  filter. 

Table  V shows  the  total  mean  error  squared  in  the  five 
states  of  the  filter.  These  numbers  were  calculated  by 
using  Equation  (67) . 


Total  Mean  Error  Squared  = j 


-2 


(67) 


t=0 


where 


= ensemble  average  of  the  state  error  over 
twenty  simulations  made  by  the  filter  at 
time  equal  to  t 
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Table 


indicates  unchanged  from  nominal  values; 


tj  = final  time  of  the  trajectory. 

The  limits  on  the  proportional  navigation  constant  were 
chosen  based  on  the  discussion  by  Stallard  (Ref  11:13-18). 
Stallard  suggested  that  for  non-optimal  guidance  the  propor- 
tional navigation  constant  should  remain  between  3 and  5. 
This  analysis  considered  6 as  a more  reasonable  upper  limit 
for  the  sensitivity  analysis.  The  range  of  the  M/S  param- 
eter was  varied  by  changing  the  mass  alone.  Since  the  two 
variables  M and  S always  appear  together  as  the  ratio 
JVS,  changing  either  will  produce  the  needed  change  in  M/S. 
The  range  of  2 slugs  to  8 slugs  was  considered  to  be  a 
large  enough  variation  to  test  the  sensitivity  of  the  filter 
to  this  parameter.  The  minimum  value  of  was  chosen  as 
0.1  seconds  in  consideration  of  the  requirement  for  the 
missile  to  filter  the  noisy  LOS  angular  rate  measurement 
(Ref  11:1).  The  maximum  value  was  chosen  at  0.8  seconds 
since  missiles  using  more  sluggish  guidance  systems  were 
considered  impractical. 

It  can  be  observed  from  Table  V that  a greater  error 
is  made  by  the  filter  at  one  extreme  of  each  parameter 
compared  to  the  other  extreme  of  the  same  parameter.  This 
sensitivity  was  not  as  noticeable  in  the  M/S  state  as  in 
the  T2  and  n channels.  The  proportional  navigation  constant 
was  varied  equally  in  both  directions  yet  the  larger  error 
committed  when  n was  set  to  3,  suggests  a higher  sensitivity 
in  this  direction.  Similar  comparisons  can  be  made  with 
respect  to  the  other  two  parameters. 
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These  observations  can  be  generalized  into  a simple 
relationship  between  the  truth  model  and  the  filter.  When 
the  parameters  of  the  filter  guidance  strip  are  equal  to 
those  in  the  truth  model,  the  dynamics  model  more  accurately 
accounts  for  the  missile's  actual  response  to  changes 
observed  in  0 . Those  changes  in  the  parameters  which  caused 
an  overall  decrease  in  the  bandwidth  of  the  guidance  strip 
in  the  truth  model,  resulted  in  larger  estimation  errors 
by  the  filter.  Conversely,  those  changes  which  increased 
the  bandwidth  caused  relatively  smaller  errors.  This 
bandwidth  comparison  can  also  be  stated  in  a lead-lag 
relationship;  more  error  is  made  in  estimation  when  the 
truth  model  lags  the  filter,  as  opposed  to  the  truth 
model  leading  the  filter.  This  relationship  can  be  observed 
in  the  zero  order  filter  developed  by  Lutter.  The  approxi- 
mation used  in  the  filter  modeled  the  guidance  strip  with 
Infinite  bandwidth.  After  this  analysis,  it  is  obvious  that 
this  was  the  primary  reason  for  the  poor  performance  of 
this  filter.  This  is  also  true  for  the  second  and  third 
order  designs  discussed  in  Chapter  III. 

The  "bandwidth  relationship"  allows  for  a heuristic 
explanation  as  to  why  the  error  is  greater  when  the  truth 
model  lags  the  filter.  Figure  26  shows  Bode  plots  of  the 
two  possible  relationships  between  the  truth  model  and 
filter  model  guidance  strips.  These  plots  were  not  calcu- 
lated but  serve  only  to  demonstrate  the  main  point  of  the 
argument.  The  relatively  wider  bandwidth  of  the  filter 


guidance  strip  model  when  the  truth  model  is  lagging,  makes 
the  filter's  model  of  the  missile  susceptible  to  more  noise 
euid  higher  frequency  information.  Since  this  high  frequency 
information  is  not  processed  by  the  actual  missile,  it  too 
is  considered  noise  because  of  its  effect  upon  the  filter 
estimate.  This  added  noise  and  high  frequency  information 
makes  the  filter  missile  appear  more  erratic  when  predicting 
forward  to  the  next  measurement.  However,  two  important 
points  lessen  the  severity  of  the  bandwidth  mismatch  when 
the  truth  model  leads  the  filter: 

1)  most  of  the  important  information  is  contained 
in  the  lower  frequencies  allowing  the  filter 
to  predict  a fairly  accurate  trajectory  and, 

2)  the  filter  is  observing  less  "noise"  in  its 
model  of  the  guidance  strip. 

The  total  error  listed  in  Table  V indicated  the  filter 
was  sensitive  to  changes  in  all  three  parameters.  This 
fact  was  also  supported  by  the  plots  of  the  filter's  error 
for  four  of  the  five  states  which  were  generated  for  the 
sensitivity  analysis.  The  complete  set  of  these  plots  are 
contained  in  Appendix  B of  this  report. 

The  results  also  indicated  that  the  filter  was  not  as 
sensitive  to  changes  in  M/S  as  it  was  to  n and  T2.  The 
filter  appears  to  be  equally  sensitive  to  n and  T2.  This 
study  did  not  consider  possible  limitations  upon  filter 
structure  during  implementation  which  would  restrict  the 


82 


mi, 


Fig,  26.  Truth  Model  and  Filter  Relationships 

number  of  filter  states.  If  such  were  the  case,  the  need 
for  estimating  n and  T2  should  be  given  higher  priority 
them  that  of  M/S. 

The  study  was,  therefore,  continued  with  the  intent  of 
estimating  all  three  parameters.  It  is  worth  noting  that, 
if  all  the  states  were  as  insensitive  to  the  parameter 
changes  as  the  range  state  is,  a nominal  value  would  have 
been  chosen  as  a result  of  this  analysis. 

Parameter  Estimation 

The  decision  to  implement  a filter  which  estimates  these 
pareuneters  would  be  based  upon  the  improved  performance 
observed  in  the  state  estimates.  The  first  order  filter  is 
the  most  likely  to  be  implemented  because  of  its  reduced 
complexity  and  general  nature.  The  filter  model  attempts  to 
approximate  the  dominant  aspects  of  an  air-to-air  missile 
without  requiring  particular  knowledge  about  the  guidance 


strip  structure.  Instead,  it  models  this  as  a first  order 
lag  which  adequately  describes  the  dominant  aspects  of  the 
guidance  strip.  The  estimation  of  n,  M/S  and  if 

estimated  correctly,  help  define  the  overall  structure  of 
the  particular  missile  being  tracked.  This  aids  the  filter 
algorithm  in  gleanir  the  important  information  from  the 
three  measurements  6^,  R,  and  R.  With  estimation  of  all 
three  parameters,  the  first  order  filter  estimates  are 
defined  belov/: 

Dynanlc  Stater 

mx  - x-component  of  missile's  inertial  velocity 

6,j,  - LOS  angle 

R - range 

R - range  rate 

- developed  lateral  acceleration  of 
the  missile 

Parameters 

n - proportional  navigation  constant 

Tj  - time  constant  of  first  order  lag 

M/S  - ratio  of  mass  to  cross  sectional  area. 

Single  Parameter  Estimation . The  first  step  in  the 
parameter  estimarion  analysis  was  designee  to  test  the 
recovery  characteristics  of  the  filter  to  initial  errors 
in  each  parameter  separately.  This  was  accomplished  by 
estimating  the  parameter  which  was  initialized  at  its  extreme, 
along  v/ith  the  five  states  of  the  missiles.  The  other  tv.’o 
parameters  were  set  at  their  correct  nominal  value.  For 
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completeness,  the  recovery  was  observed  from  both  extremes 
of  each  parameter.  The  extremes  for  the  navigation  constant 
were  kept  the  same  as  established  in  the  sensitivity  analy- 
sis, ranging  from  3.0  to  6.0.  The  M/S  ratio  was  varied  from 
2 2 

15.0  slugs/ft  to  45.0  slugs/ft  with  the  true  value  at  29.197 

2 

slugs/ft  . This  corresponds  to  a change  in  the  mass  from 

2.0  slugs  to  6.0  slugs.  The  parameter  was  varied  from  0.3 
to  1.5  seconds.  The  range  for  was  based  upon  the  .85 
seconds  value,  which  resulted  in  the  best  estimation  when 
matching  the  first  order  lag  to  the  complete  guidance  strip 
(Reference  Chapter  III)  . Figures  27  through  32  show  the 
recovery  characteristics  of  the  three  parameter  estimates  from 
poor  initial  values.  Two  important  observations  can  be  made 
concerning  these  plots ; 

1)  All  three  parameters  can  be  estimated  separately. 

2)  All  three  parameters  require  an  acquisition  time 
of  approximately  2 seconds  (the  acquisition  time 
is  defined  as  the  time  necessary  for  the  mean  to 
approach  and  remain  within  one  sigma  of  its 
apparent  steady  state  value) . 

Also  observed  during  this  procedure  was  the  improved  estima- 
tion of  the  five  states  after  the  acquisition  time  (i.e., 
once  adequate  parameter  estimation  had  been  obtained) . The 
best  example  of  this  is  velocity  plot  shown  in  Figure 

33.  This  plot  was  made  while  estimating  t^,  which  was 
initialized  at  0.3  seconds.  The  plot  can  be  compared  to 
Figure  34  which  shows  the  velocity  which  was  generated 
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Fig.  33.  X-VELOCITY  FIRST  ORDER  MISSILE 
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with  perfect  knowledge  of  all  three  parameters. 

Simultaneous  Estimation  of  n and  After  having 

satisfactorily  estimated  each  parameter  separately,  there 
was  still  some  question  as  to  whether  good  state  estimation 
was  possible  when  estimating  more  than  one  parameter  at  a 
time.  This  question  was  partially  answered  by  simultaneous 
estimation  of  n and  t^.  This  combination  of  two  parameters 
was  considered  to  be  the  most  difficult  to  estimate  together 
since  n and  both  appeared  only  in  the  guidance  strip  of 
the  filter  model.  Also,  the  variation  of  each  of  these  param- 
eters strongly  influenced  the  bandwidth  of  the  guidance 
strip.  Thus,  simultaneous  estimation  of  n and  addressed 
the  question  of  whether  the  filter  could  distinguish  between 
the  changes  needed  to  converge  to  a proper  estimate  and 
those  changes  possibly  made  to  compensate  for  the  effect 
of  the  other  parameter. 

The  answer  to  this  question  can  best  be  given  by 
showing  the  estimation  plot  (Figure  35)  . The  acquisi- 

tion time  for  the  navigation  constant  and  the  time  constant 
for  the  first  order  lag  (shown  in  Figure  36  and  37) 
remained  at  2 seconds.  After  this  time,  the  greatest 

error  observed  in  the  channel  is  40  ft/sec.  This 

mx 

compares  favorably  with  the  result  of  35  ft/sec  reported  in 
Chapter  III,  when  the  parameters  are  assumed  perfectly 
known. 

Up  until  this  point  in  the  study,  all  simulations  had 
been  accomplished  using  the  high-g  scenario.  There  was 
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doxibt  as  to  whether  the  filter  would  yield  accurate  state 
estimation  when  also  estimating  the  parameters  for  the  low-g 
trajectory.  The  doubt  seemed  justified  for  the  following 
reason.  Although  the  parameters  were  not  functions  of  the 
missile  states,  their  estimation  was  dependent  upon  the 
excitation  of  these  states.  The  wider  the  range  of  varia- 
tions observed  in  the  missile  states,  the  more  information 
the  filter  was  capable  of  gleaning  from  the  measurements 
for  estimating  the  parameters.  The  low-g  trajectory  was 
benign  in  comparison  to  the  high-g  scenario  and  offered 
minimal  state  variations. 

Figure  38  confirms  the  suspected  parameter  estimation 
problem  for  the  low-g  trajectory.  The  best  example  of 
this  problem  was  the  parameter.  Without  adequate 
parameter  estimation,  it  seemed  less  likely  that  enhanced 
state  estimation  was  possible;  in  fact,  performance  might 
even  be  degraded  from  that  of  a filter  with  a less  refined 
acceleration  model.  Yet,  Figure  39  shows  this  statement 
to  be  false,  since  the  error  was  less  than  30  ft/sec 
for  the  entire  scenario.  This  result  indicates  that, 
although  beneficial  for  the  high-g  scenario,  the  refined 
acceleration  model  was  not  advantageous  for  the  low-g 
trajectory.  This  also  implies  that  a Gauss-Markov  accelera- 
tion model  could  possibly  be  used  for  low-g  scenarios. 
Nevertheless , the  refined  acceleration  model  did  not  intro- 
duce severe  estimation  problems  for  the  low-g  case,  and 
is  therefore  considered  generally  applicable  over  a wide 
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range  of  encounter  trajectories. 


Siroultaneous  Estimation  of  All  Three  Parameters.  The 
success  of  estimating  n and  together  was  followed  by 
the  logical  decision  to  estimate  all  three  parameters 
simultaneously.  The  high-g  scenario  was  used  to  investi- 
gate the  parameter  estimation  better.  Figure  40  illustrates 
the  important  points  of  this  estimation.  The  other  plots 
are  contained  in  Appendix  B. 

The  estimate  of  each  parameter  was  not  appreciably 
degraded  with  the  increased  complexity  of  estimating  all 
three.  The  acquisition  time  of  about  2 secs  can  be  observed 
in  the  error  plot.  The  maximum  error  after  this  time 
is  less  than  50  ft/sec.  The  increasing  error  in  the  first 
2 secs  of  the  plot  is  attributed  to  the  initial  error  in 
the  parameters . 

Summary 

A sensitivity  analysis  was  performed  which  indicated 
that  the  three  key  parameters,  n,  M/S  and  i could  be 
estimated  along  with  the  five  states  of  the  missile.  The 
large  errors  committed  by  the  fourth  order  filter,  during 
the  sensitivity  analysis,  strongly  suggested  that  choosing 
a nominal  value  for  any  parameter  would  significantly 
degrade  performance.  This  is  especially  true  for  the  first 
order  filter  which  did  not  include  as  refined  a model  of  the 
guidance  strip  as  that  of  the  fourth  order  filter.  Because 
the  tradeoff  was  considered  so  great,  no  suggested  nominal 
parameter  values  are  included  in  this  study.  A 
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systematic  estimation  analysis  was  then  performed  which 
strongly  confirmed  that  simultaneous  estimation  of  all  three 
parameters  and  states  was  possible.  In  order  to  diminish 
the  transient  effects  of  the  states  and,  to  observe  the 
effects  of  parameter  estimation  on  the  states  better, 
the  dynamic  states  were  initialized  with  zero  error.  This, 
of  course,  is  unrealizable  in  a real  world  situation.  The 
final  test  of  the  filter  was  to  subject  it  to  estimation  of 
all  five  states  and  all  three  parameters.  The  states  were 
initialized  with  some  reasonable  error  (10%  off  the  true 
value)  and  the  parameters  were  initialized  at  one  of  their 
extremes.  The  tuning  parameters  used  when  initializing 
the  states  with  zero  error  were  found  inappropriate  for 
this  case.  This  also  implies  the  possible  need  for  adaptive 
self  tuning  in  a filter  considered  for  implementation 
since  the  performance  presented  in  Chapter  V would  not  have 
been  achieved  without  ratuning. 
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V.  Results  and  Conclusions 

Introduction 

A major  purpose  of  this  thesis  has  been  to  design  a 
reduced  order  filter  which  could  estimate  the  dyncunic  states 
euid  uncertain  parameters  of  an  air-to-air  missile  model. 

The  first  order  missile  model  filter  developed  in  the  last 
chapter  has  great  potential  for  this  application.  However, 
the  dedication  to  this  effort  of  estimating  all  three 
uncertain  parameters  in  this  model's  structure  may  have,  up 
to  this  point,  overshadowed  the  prospects  of  only  estimating 
one  or  two  of  these  parameters.  Although  Chapter  IV 
emphatically  states  that  the  performance  degradation  is  great 
for  a "mismatched"  nominal  parameter  value,  it  is  possible 
that  the  composition  of  a particular  problem  may  either  have 
considerably  smaller  ranges  for  the  values  of  the  nominals, 
or  be  constrained  in  some  way  from  estimating  all  three 
parcuneters.  For  this  reason,  this  chapter  begins  with  a 
discussion  on  estimating  fewer  than  three  uncertain  param- 
eters. 

The  final  evaluation  of  the  first  order  missile  filter 
is  then  presented.  This  evaluation  attempted  to  estimate 
simultaneously  all  three  parameters  and  all  five  states.  The 
final  evaluation  is  distinguished  from  previous  cases  by  the 
fact  that  the  states  as  well  as  the  parameters  are  initialized 
with  some  error.  The  results,  presented  here,  continued  to 
demonstrate  the  versatile  nature  of  the  first  order  missile 
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filter.  The  significant  findings  of  this  study,  as  well  as 
suggestions  for  future  work,  are  presented  following  these 
results. 


Table  VI 

Total  Error  in  Each  State 
During  Single  Parameter  Estimation 


^\^tate 

Pa- 

raiaete^v 

Estimated^ 

vL  (10^) 

UK 

e do"'*) 

R (10^) 

R (10^) 

\ (10^) 

none 

0.99 

3.86 

1.82 

0.478 

3.70 

JVS 

1.39 

5.21 

1.85 

1.870 

5.93 

n 

2.84 

15.20 

1.85 

0.514 

7.23 

^f 

3.45 

15.50 

1.94 

1.370 

1 

11.10 

1 

Estimating  Fewer  Than  Three  Uncertain  Parametv  rs  i 

Table  VI  presents  the  total  mean  error  squared  in  the 
dynamic  states  of  the  filter  model  (calculated  according  to  j 

the  method  discussed  in  Chapter  IV)  for  the  estimation  of 
single  parameters.  The  initial  entry  shows  the  error  gen- 
erated with  perfect  knowledge  of  these  parameters.  The 
parameter  that  was  estimated  is  listed  in  the  left  hand 
column.  This  parameter  was  initialized  with  some  error  while 
the  other  two  parameters  were  set  at  their  correct  values. 

The  total  errors  when  estimating  the  individual 
parameters  show  that  the  filter  is  more  sensitive  to 
initial  errors  in  the  and  n parameters  than  in  M/S. 

This  observation  is  consistent  with  the  sensitivity 
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analysis  performed  on  the  fourth  order  missile  filter  which 
indicated  that  the  filter  was  more  sensitive  to  T2  and  n 
compared  to  M/S'.  Table  VI  delineates  more  fully  the 
priority  for  including  the  parameters  in  the  first  order 
missile  model.  Therefore,  if  confronted  with  the  requirement 
to  choose  a nominal  for  one  or  more  of  these  parameters,  it 
is  suggested  that  the  M/S  parameter  be  given  first  consider- 
ation. The  table  indicates  that  the  n parameter  should 
be  considered  next  with  chosen  last.  However,  because  of 
the  small  difference  between  the  total  errors  listed  for  n 
and  T^,  a particular  problem  may  justify  the  estimation  of 
n as  opposed  to  x^.  The  choice  of  nominal  values  is  very 
much  problem  dependent,  and  because  of  the  tradeoff  in 
filter  performance  when  using  nominals,  no  values  have  been 
recommended  by  this  study. 

Final  Results ; Estimating  All  Uncertain  Parameters 

Figures  41  through  48  exhibit  the  performance  of  the 
five  states  and  three  parameters  for  the  final  evaluation. 
Table  VII  presents  the  key  results  along  with  information 
from  the  benchmark  design  and  fourth  order  filter  for  com- 
parison. The  initial  errors  for  all  missile  states  were 
considered  within  the  performance  ability  of  a typical 
fighter  aircraft  radar.  The  V^  error,  which  v^as  initial - 
ized  with  an  error  of  225  ft/sec,  indicates  a recovery 
time  of  one  second.  The  recovery  time  is  defined  to  be  the 
time  necessary  for  the  initial  error  induced  transient  to 
decay  to  the  performance  produced  in  a filter  initialized  with 
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zero  dynamic  state  error.  The  maximum  error  after  this 
recovery  time  was  about  40  ft/sec.  This  can  be  compared  to 
the  35  ft/sec  error  of  the  first  order  missile  filter  with 
perfect  knowledge  of  all  three  parameters. 

The  angle  error,  which  was  initialized  with  an  error  of 
.02  radians  (1.1  degrees),  indicates  a recovery  time  of  1.7 
seconds.  This  initial  error  can  be  considered  as  a cross 
rauige  error  of  approximately  180  ft  at  the  start  of  the 
intercept.  The  maximum  error  of  this  state  never  exceeds 
2 milliradians  after  the  engagement  has  progressed  past  the 
recovery  time. 

The  range  and  range  rate  errors  exhibit  very  fast  recovery 
times  of  .5  seconds.  The  errors  committed  after  this  time 
were  extremely  small.  These  two  results  were  due  to  the  fact 
that  both  of  these  states  are  receiving  information  from 
accurate  direct  measurements. 

The  lateral  acceleration  state  error  shows  the  same  1.7 
second  recovery  time  shown  in  the  angle  state.  It  was 
initialized  with  an  error  of  15  g's.  Although  this  may  be 
considered  an  unreasonably  large  error  for  most  cases,  it 
could  possibly  arise  while  tracking  a missile  soon  after 
launch.  In  this  case  the  missile  may  enter  a high-g 
maneuver  in  an  attempt  to  null  initial  launch  errors.  It 
is  worth  noting  that,  at  five  seconds,  the  missile  is 
achieving  a lateral  acceleration  in  excess  of  15  g's. 

However,  the  error  in  the  filter  estimate  at  this  time  is 
only  1 g.  Again,  this  is  evidence  of  very  good  state 
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estimation. 

The  filter  does  not  acquire  an  accurate  estimate  of  the 

proportional  navigation  constant  until  approximately  2.5 

seconds  after  initialization.  Subsequent  to  this  time,  the 

filter  mean  error  remains  less  than  5%  of  true  value  (4.5). 

The  parameter  estimate  does  not  converge  to  zero  error 

as  well  as  n does,  but  shows  continual  improvement  over  the 

entire  scenario.  After  2.5  seconds,  it  was  estimated  well 

enough  to  aid  in  accurate  state  estimation.  The  error 

in  the  M/S  ratio  reaches  steady  state  the  soonest  of  all 

three  parameters.  It  converges  quickly  (approximately  one 

2 

second)  but  remains  biased  by  a positive  2 slugs/ft  . 

This  is  only  6%  of  the  true  value  (29.197  slugs/ft).  These 
results  confirm  that  the  first  order  missile  filter  is 
capable  of  estimating  both  the  states  and  parameters  of  an 
attacking  air-to-air  missile. 

Contributions 

The  contributions  of  this  study  include  the  following. 
Given  the  three  independent,  radar-supplied  measurements  of 
LOS  angle,  range,  and  range  rate  of  an  air-to-air  missile, 
an  eight-state  Extended  Kalman  Filter  (modeling  five  dynamic 
states  and  three  key  parameters  of  the  missile)  was 
designed  such  that: 

la)  Simultaneous  estimation  of  the  states  and 
parameters,  with  reasonably  accurate 
results,  was  possible  for  a high-g  scenario. 

lb)  The  filter  was  considered  "robust"  since,  when 
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tracking  a low-g  target  the  filter  was  capable 
of  achieving  good  state  estimation  even  with  the 
parameter  estimation  not  being  as  accurate  as  for 
the  high-g  scenario 

2)  Modeling  efforts  for  this  type  of  filter  must 
consider  the  bandwidth  of  the  guidance  strip. 

This  fact  can  be  emphasized  by  recalling  that 
even  including  an  exact  model  of  the  missile's 
autopilot  did  not  signficantly  aid  estimation. 

3)  The  filter  could  recover  from  large  initial 
errors  in  the  dynamic  states  and  parameters  of 
the  missile. 

4)  Implementation  of  this  filter  would  require  an 
acquisition  time  to  lock  on  to  adequate 
parameter  estimates.  The  best  example  of  this 
was  the  2.5  second  acquisition  time  of  propor- 
tional navigation  constant. 

Future  Research 

The  following  suggestions  are  listed  as  possible  areas  for 
future  research: 

1)  The  analysis  considered  only  non-thrusting 
missiles.  It  is  suggested  that  this  restriction 
could  be  removed  by  using  a random  walk  to  model 
the  acceleration  due  to  thrust. 

2)  Because  all  filters  were  optimally  tuned  for  a 
given  trajectory,  the  robustness  observed  over 
that  one  trajectory  may  not  be  observed  for 
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other  scenarios  without  retuning  the  filter. 

This  suggests  that  an  analysis  be  peformed  which 
more  completely  investigates  the  "robustness" 
of  the  filter,  possibly  leading  to  the  use  of 
adaptive  tuning  techniques  in  the  filter  design. 

3)  Parameter  estimation  was  not  needed  for  the  low-g 
trajectory  to  achieve  accurate  state  estimates. 
The  initial  trajectory  of  an  air-to-air  missile 
may  in  some  cases  be  this  benign  since  the  LOS 
rate  is  small  for  extended  ranges.  This  suggests 
the  possibility  of  an  adaptive  filter  which 
initially  uses  a random  walk  or  first  order 
Markov  model  for  acceleration  and  then  switches 
to  the  more  complete  model  presented  in  this 
report. 

4)  The  analysis  should  be  expanded  to  consider 

- - - A 

three  dimensional  motion  of  the  missile  and 
the  ownship. 

5)  An  analysis  could  be  performed  which  investi- 
gates the  degradation  of  filter  performance  by 
increasing  the  measurement  update  time  interval. 
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Appendix  A 


This  appendix  contains  the  changes  to  the  software 
package  designed  by  Lutter  as  reported  in  Reference  5.  There 
were  five  changes  to  the  original  software  which  resulted 
in  software  used  for  this  study.  The  changes  were: 

1)  The  complexity  of  the  integration  package  of  the 
filter  was  reduced  to  a first  order  Euler. 

2)  The  initialization  of  the  matrix  subtraction  sub- 
routine was  changed. 

3)  The  modeling  of  the  seeker  noise  in  the  truth  model 
was  changed  to  include  both  correlated  and  uncor- 
related noises. 

4)  The  original  noise  strengths  were  scaled  (multi- 
plied) by  a ,it  factor. 

5)  The  updating  of  the  lineraized  F matrix  was  accom- 
plished at  each  integration  step. 

Also  included  is  a discussion  on  the  assumption  made  by 
Lutter  concerning  the  linearization  of  the  F matrix 
tRef  5:65) . 

Euler  Integration  Package 

The  reduced  complexity  of  the  integration  package  of  i 

i 

the  filter  was  accomplished  by  replacing  a fourth  order  ' 

1 

Runge-Kutta  package  with  a first  order  Euler  integration 
package.  The  Runge-Kutta  package  was  constructed  with  a 
backward  predictor  corrector  which  automatically  decreased 
the  integration  step  size  if  the  backward  prediction  was 
not  within  a specified  tolerance  of  the  starting  value. 
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This  one  change  reduced  the  execution  time  of  the  simulation 
to  one-fifth  of  its  original  time. 

Even  though  a simple  integration  proceedure,  the  Euler 
package  did  not  degrade  the  simulation  accuracy  substantially 
A sample  case  was  run  using  the  Runge-Kutta  package  and 
the  Euler  package  for  the  same  data  set.  The  plots  of  the 
case  using  the  Euler  integration  pacakage  were  indistin- 
guishable from  the  plots  generated  using  the  Runge-Kutta 
package.  These  results  are  shown  in  Appendix  B.  The  Euler 
package  was  also  used  because  it  more  closely  approximated 
the  integration  packages  available  for  on-line  implementa- 
tion purposes. 

Chanae  Matrix  Subtraction 


The  purpose  of  the  matrix  subtraction  subroutine  was 
to  perform  the  subtraction  operation 


C = A + B 


where  A,  B,  and  C were  all  matrices  of  the  same  size. 

The  original  routine  was  set  up  so  that  the  C matrix  was 
initialized  at  zero.  This  was  unnecessary  because  C was 
not  used  in  the  actual  calculation,  but  just  as  a location 
to  store  the  results.  The  original  routine  also  led  to 
problems  when  passing  information  betv;een  the  subroutine 
and  the  main  program.  The  same  argument  was  used  for  the 
A matrix,  passed  into  the  subroutine,  and  the  resultant 
matrix  C,  passed  back  to  the  main  program,  within  the 
subtraction  subroutine.  This  resulted  in  setting  all 


values  in  the  matrix  A to  zero.  Thus 

C = 0 -B  = -B  (59) 

To  avoid  the  problem  and  to  simplify  the  subroutine  the 
zero  initialization  of  C was  removed. 

Modeling  Seeker  Noise 

In  his  software,  Lutter  inappropriately  modeled  the  seeker 

noise  so  as  to  include  only  a time-correlated  noise.  As 

stated  in  Chapter  II,  to  model  this  noise  correctly,  both 

the  correlated  and  uncorrelated  components  were  included. 

The  At  Factor 

When  developing  an  equivalent  discrete  system  of  equa- 
tions for  the  Extended  Kalman  Filter,  noise  strength  must 
be  given  special  consideration.  The  noise  sigmas,  given  in 
Chapter  II,  are  related  to  continuous  noise  sources.  To 
model  this  continuous  noise  correctly  by  discrete  equations, 
the  noise  strengths  should  be  multiplied  by  a factor  of  At, 
the  integration  time.  This  At  factor  was  included  in  this 
study . 

To  explain  the  reason  for  including  the  At  factor,  the 
equivalent  continuous  and  discrete  time  noise  relationships 
should  be  compared.  Consider  a linear  stochastic  differen- 
tial equation  of  the  form; 

x(t)  = F(t)x(t)  + G(t)w(t)  (70) 

where 

w(t)  = a zero  mean  white  Gaussian  noise  with 
E{w(t)w(t+T)’^}  =^(t)6(t) 
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An  equivalent  discrete  time  system  equation  would  be 

t . , 

^^^i+1^  = i(t^+i't.)x(t^)  + /^^■^■^i(t^^^,t^)G(T)w(T)dT  (71) 


= i(ti+i.t.)x(t.)  w^(t.) 


where 


E{Wj(tj^)W^(t^)’')  = /t^'^'’i(ti^,,T)G(_t)Q^(T)G(T)’'»  (tj  , ,T)’'dT 


= ea(ti) 


E{w,  (t . ( . )'^}  = 0 for  t.  f t. 

-di-dj  — 


When  the  integration  step  size  is  small  relative  to 
system  time  and  constants,  the  following  approximation 


can  be  made 

^i+1 

/ ^ ■ £G0  cVdT  = G(t.)0  (t.)G(t,)'^At  (7; 

2 

where  all  terms  of  order  At  or  higher  have  been  neglected. 


Therefore 


Q^(t.)  = {G(t^)Q^(t^)G(t^)MAt 


and  for  G(tj^)  = I (the  identify  matrix) 


Q^(ti)  = Q^Ct.)At 


This  argument  is  also  directly  applicable  to  nonlinear 
stochastic  differential  equations  in  which  the  driving 
noise  enters  linearly,  and  so  it  was  used  in  this  study. 


Update  of  the  F Matrix 

Because  of  the  complicated  and  time  consuming  nature  of 
the  Runge-Kutta  integration  package, and  the  use  of  a very 
short  measurement  update  time  of  .02  seconds.  Latter  made 
the  decision  to  update  the  linearized  F matrix  of  the 
Extended  Kalman  Filter  only  at  measurement  update  times. 

The  integration  package  in  this  study  was  less  complicated 
and  therefore  reduced  the  execution  time.  It  was  reasonable 
to  update  the  F matrix  at  each  integration  step  (.01  seconds). 
This  did  not  increase  the  execution  time  any  appreciable 
eunount  and  was  a more  correct  implementation  of  the 
Extended  Kalman  Filter  equations. 

Linearization  of  f-vector 

Lutter  assumed  that  Y / and  V , v/hich  are  all  func- 

T'  'm'  m 

tions  of  the  filter’s  state  variables,  remained  constant 
over  the  propagation  period  of  .02  seconds.  One  example 
of  how  this  assumption  affects  the  linearization  of  the  f- 
vector  can  be  shown  by  the  development  of  the  partial 

0 

derivative  of  ^ for  an  element  which  is  a function  of  9,^ 


3f  • {6_(x.  f . . • f • • • ^ri  ^ 

1 T 1 n 1 n 

3x . 


(76) 


3f  36  3f. 

1 + — £ 

36  3Xj  3Xj 


6^*  constant 


(77) 


The  assumption  considers  the  first  term  of  Equation  (77) 
to  be  zero.  This  assumption  was  verified  by  computer  simula- 
tion of  the  high-g  scenario  for  the  fourth  order  missile 
filter.  A complete  linearization  of  the  f-vector  which  did 
not  assume  these  terms  to  be  zero  was  used  in  this  filter. 

The  plots  for  all  states  are  labeled  "comp"  and  can  be  found 
in  Appendix  B.  Comparison  of  these  plots  with  the  fourth 
order  set  labeled  "A/P  EQ  0",  whose  filter  F-matrix  included 
these  assumptions,  confirms  the  validity  of  the  assumption 
made  by  Lutter. 
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